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LXXVIIT. Calibration of Proportional Counters by the Excitation of 
Flourescence Radiation with Radioactive Sources. 


By GoM. Insou; 
Department of Natural Philosophy, Glasgow University *. 


[Received July 3, 1950]: 


SUMMARY. 


The calibration of proportional counters by means of a very few suitably 
chosen radioactive sources which are used to excite the fluorescence 
radiation of various elements is described. The efficiency of excitation is 
high and the method is of wide application. The accurate identification 
of isomeric transitions using proportional tubes is demonstrated. 


INTRODUCTION. 


THe Investigations carried out in this laboratory using proportional 
counters (Phil. Mag. 1949) have necessitated detailed examination of various 
calibration techniques. Information on this subject has already been 


Big: 1. 
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published by Rothwell and West (1950). However, during the course of 
our work considerable progress has been made in overcoming the various 
problems involved, and it is the aim of this paper to correlate some of the 
results obtained. The chief method described here is original and it is 
described in reasonable detail since it is of wide application. 

In the early investigations the beam from an ordinary A.C.-operated 
X-ray machine, adjusted to operate at a small output intensity and which 
emitted white radiation was allowed to fall on a thin foil and the scattered 
radiation was detected in the counter at an angle to the beam direction 
(fig. 1(a)), (Curran, Angus and Cockroft,1949). This radiation entering the 


* Communicated by the Author. 
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tube consisted mainly of the fluorescence X-rays of the scattering element. 
However, this method is not ideal. The X-ray machine is cumbersome, 
may give rise to spurious pulses (due to electrical interference), and calibra- 
tion is accompanied by an increase in the “ background ” of the counter 
due to the detection of the white radiation. In addition the tube in ageing 
fails to self-rectify adequately and hence the pulses appear grouped on part 
of alternate half-cycles of the mains supply. For these reasons, new 
calibration techniques were explored and normally preferred. Pile- 
activated sources, which decayed by K-capture, or which showed nuclear 
isomersion or strong internal conversion of their y-radiation, were found 
suitable as agents for the production of beams of homogeneous X-radiation 
or electrons. Sometimes the sources were used directly, but two methods 
of greatly entending their usefulness were established. The sources which 
have been most commonly applied in the course of several studies are shown 
in Table I. 


TABLE I. 
Teoteuen? sHalbifems “trdneitign co PP Os CNCTEY es eee 
of X-ray radiation 
Tn 186 50 days isomeric 24-1 KeV. 1-98 MeV. 8+y’s 
Snis 105 days K-capture 24-1 KeV. y's and e- 
Gere 11 days K-capture 0-2 Ke. negligible 
Fe:8 4 years K-capture 5-9 KeV. none if pure. 


ELIMINATION OF THE Hrrects or Harp Beta AND GAMMA RADIATION. 


The sources should have a reasonable half life for prolonged use. They 
should preferably be thin, particularly when associated with hard y-radi- 
ation. Experience has shown that isotopes which have long period 
B-activity may be comfortably handled. The f-particles are absorbed 
by suitable thicknesses of polythene which has good transparency for 
X-radiation. For example, in the investigation of In!4 (see below) the 
source was covered by a layer of polythene, 6mm. thick. This completely 
absorbed the hard f-radiation (of energy 1-98 MeV.) and only reduced the 
X-radiation (energy 24:1 KeV.) by a relatively small fraction (~1/4). 
When the required soft calibrating X-radiation is associated with 
y-radiation (of energy say ~100 KeV.) some shielding is useful in reducing 
the undesirable y-ray flux through.the counter to a miniumm (fig. 1(b)). 
However, this method can only be applied with advantage to relatively soft 
y-rays. In the case of still harder y-radiation (~-25 MeV. upwards) the 
source is most conveniently examined when placed really close to the 
counter window. ‘This simple geometry results in the best solid angle for 
soft radiation passing into the counter and makes use of the higher efficiency 
of the counter for such radiation. Usually it is possible to secure a ratio 
of at least 10 : 1 in efficiency of detection in favour of the soft X-radiation. 
To obtain optimum results when examining the source in this way, care 
must be taken to eliminate the “‘ end-effect ” of the tube, since a small 
fraction of the radiation detected penetrates to the region of variable gas 
gain. 
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Propuction or Mopirtep Raprarion. 


Calibrating agents which gave homogeneous X-radiation covering the 
whole range from say 3 KeV. to 50 KeV. in small energy increments were 
found extremely useful. Since this coverage cannot be accomplished 
readily by use of radioactive sources alone new methods were devised. 
These involved the excitation of the desired radiation by homogeneous 
X-rays of higher energy originating from a radioactive source. This 
process is one of high efficiency and it was achieved in practice either by 
allowing the initial radiation to fall on a thin sheet of the appropriate 
element and detecting the excited radiation emitted from the front face 
(reflection (fig. 1(c)), or by placing a suitable thickness of the element 
between the source and the window (transmission). In the former case, 
the fluorescence radiation greatly predominates in intensity within the 
counter and its homogeneity is nearly independent of the thickness of the 
reflecting sheet. In the second case (transmission) a suitable thickness of 
the sheet can readily be calculated. For example, consider a source emitting 
homogeneous X-rays of wavelength A, exciting a radiation of wavelength 
A, in a sheet of the chosen element placed between the source and the 
counter window. Let this sheet have mass absorption coefficients j., 
and jz. corresponding to A, and A,. Let it be of thickness / and density p. 

It can be shown that the number of fluorescence X-ray quanta entering 
the counter (fig. 1(d2)) when I, quanta fall on the sheet is given by 


N=K,I,[e-™*—e-™]/(K,—K,), 
where 
Ko, pe and. Kap, pl. 


The number of unmodified X-rays transmitted through the sheet to the 
counter is 


and the ratio of these is 
R=N/T}K fee! — 1) (Ki). 2 oe a 0) 


Examination of equation (1) shows that the ratio R can be made as large 
as required (>100 say) by using sufficient thickness J, provided that 
y<p4,, a condition which is always satisfied in the usual experimental 
arrangements. (Furthermore this argument is strengthened by the fact that 
the modified radiation is usually the more efficiently detected). However, 
in practice large / values may necessitate strong sources and in some cases 

it may be difficult to secure adequate intensity. On the other hand, if R~1 
is regarded as suitable (peaks due to modified and unmodified radiations of 
equal intensity) the method can be used with very weak sources. (see 
fig. 2). . . a . 

Regarding the reflection method (fig. 1(¢)) it 1s sufficient to note that it 
can be readily applied in almost every case, but it may be inferior to the 
transmission method when very weak sources are employed. 


BO? 
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These methods were applied to excite the fluorescence radiations of a 
number of elements (see Table II.). In this Table E, is the approximate 
energy of the initial X-radiation falling on the element and E,, is the 
approximate energy of the modified radiation. 

With the simple geometry demonstrated in fig. I(c) it was possible to 
detect a) modified radiation ~5 per cent of the total radiation emanating 


Fig. 2. 


15 7 Kev. 


Number of pulses per energy interval. 


Energy in arbitrary units. 


from the source. Hence a source strength of <I curie total activity 
proved sufficient. Indium and Geranium sources were found adequate 
for excitation of the flourescence radiations from less than 5 KeV. upwards, 
which considerably increased the number of homogeneous X-ray 
calibrating agents available within this energy range. 


Tasie II. 
ewe 
Source eee KE, in KeV. E,, eid ie Histogram 
ines Zr 24-1 15-7 Tramsnission fig. 2 
Inits Cu 24-1 8-05 Reflection fig. 3(a) 
Ges Ni 9:2 75 Reflection fig. 3(b) 
Ge Fe 9:2 6-4 Reflection fig. 3(c) 


New calibrating agents are at present being examined for use with 
proportional counters in magnetic fields. Curran, Cockroft and Insch 
(1950) have recently used proportional tubes to measure particle energies 
>1MeV. These examinations demand calibration radiations of several 
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hundred KeV. energy. Sources which decay by y-ray transitions (not 
associated with 6-aetivity) and which show strong internal conversion, are 
being studied (e.g. Sn'!8; y-ray energy 392 KeV. and ~70 per cent 
internal conversion) Such sources are introduced into the counter. 


IDENTIFICATION OF ISOMERIC TRANSITIONS. 


The primary phenomena externally observable in an isomeric transition 
are the emission of y-radiation and usually conversion electrons. Accom- 
panying the latter is a K X-radiation which is detected with high efficiency 
and accuracy by a proportional counter. Since this X-radiation is 


Fig. 3. 
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frequently difficult to identify with certainty, the method described here is 
an important aid. The isomeric transition in In‘! was studied and the Ke 
fluorescence X-ray histogram is shown in fig. 4, curve A. At the same time, 
the Kx fluorescence radiations of cadmium (energy 23-1 KeV.) Indium 
(24-1 KeV.) and tin (25-2 KeV.) were excited and compared with the radi- 
ation arising in the nuclear process. The three histograms are shown in 
fig. 4, curves B, Cand D. This demonstrates at once that the known y-ray 
energy 192 KeV., which is strongly converted, is emitted prior to any 
transition in which Z changes. The importance of this method for the 
accurate identification of isomeric processes is obvious. 

Tn the course of this work argon was used as the absorbing gas. For the 
energies investigated it is superior to the gases of greater atomic number, 
such as Xenon and Krypton, which have a relatively small Auger effect, 
causing homogeneous X-radiation to yield two peaks of comparable 
intensity, which can lead to uncertainties. One peak corresponds to the 
case when the argon X-ray escapes from its own atom and is not detected. 
The Auger effect in argon is high and an independent yalue for this effect 
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was obtained during the work described above. This was accomplished 


by studying an X-ray histogram obtained when the fluorescence K- 


X-rays of iron were excited by a germanium source. ‘Two peaks were 
evident, one due to the fluorescence Ka X-radiation of iron of energy 
6-4 KeV., and one to the photoelectrons (energy 3-4 KeV.) from the K level 
of argon when the X-ray escaped and did not form an Auger electron. 


Fig. 4. 
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The effect on the smaller energy peak, due to 3 KeV. X-radiation caught 
in the gas was calculated and the “ background ” spectrum of the counter 
was subtracted from the histogram. Using this the Auger effect in argon 
was found to be 94 per cent in reasonable agreement with previous work 
(Compton and Allison 1936). 

Throughout this work a counter 54” internal diameter was used with a 
filling of 60 em. Hg. of argon and 15 cm. Hg. of methane. 


REFERENCES. 
RotrHweE ., P., and Wzst, D., 1950, Proc. Phys. Soc., 63, 541. 
Curran, 8. C., Anaus, J., and Cockrort, A. L., 1949, Phil. Mag., 40, 36. 
Curran, S. C., Cookrort, A. L., and Iyson, G. M., in the press. 
Compton, A, H., and Attison, 8. K., 1936, X-rays (Macmillan). 


[86a] 


UX XEXS Transit-time Phenomena in Electron Streams.—Il. 


By D. K. C. MacDonatp, 
Clarendon Laboratory, Parks Road, Oxford *. 


[Received July 13, 1950.] 


ABSTRACT. 

A density modulated electron stream will suffer amplitude reduction 
and phase variation with transit time due to the statistical velocity 
distribution of the electrons. The problem is considered in relation to 
that of noise generation in a beam due to velocity fluctuation; the 
analysis developed is then used in examining the variation of space-charge 
reduction factor (I) under particular laws of potential distribution 
which present themselves in physical problems. 


$1. INTRODUCTION. 

IN a previous paper in this Journal + (MacDonald 1949 a) we analysed the 
problem of an electron beam initially accelerated to a relatively high 
potential and initially ordered in density (I’=0), which becomes progres- 
sively disordered (J 1) due to the random distribution of velocities 
present in the electron beam. We wish now to consider the related 
problem of an electron stream with thermal velocities which has a given 
density-modulated sinusoidal ** signal ”’ initially imposed ;. the amplitude 
and phase of this signal will then suffer alteration with transit-time due 
to the thermal energy spread. We shall examine the significance of this 
relationship and then proceed to analyse the variation of such a signal 
and the related J under certain potential laws of physical interest. 
The theoretical predictions will be compared with such experimental 
evidence as is available and the significance of the results discussed. 


§2. REDUCTION OF SIGNAL AMPLITUDE. 

Since (I) was published, both van der Ziel (1949) and MacDonald 
(1949 b) have independently considered the variation du: ing drift time’ 7 of 
a sinusoidal density-variation in an electron stream of thermal energy 
distribution 

p(E)dE=1/k0 . exp(—H/k0)dE, (E>0) l 
=() ; (EK <0) J 
(k : Boltzmann’s constant. 0: Absolute temperature) 
which has been initially accelerated very rapidly to energy ®. 
If the beam-density before “‘ drifting ” (70) is given by 
Nqy=No+ Apo sin wi, 


(1) 


* Communicated by the Author. 
+ Referred to as (I) here, 
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then one shows readily by considering the loss and gain of electrons in a 
short interval of time that 
2D Gt = ( 20x 


Ng) aa : ex ) . (Ny +Ap sin wt cos wa-+ Ay sin wx cos wt) dx 


kOe 
=N,+A, sin ot | a} +o cos of | th, a ean utas 
ké 
where a=wt (za): 


The unvarying (“ d.c. ’’) component, NV, remains unchanged, as we should 
expect, but the oscillatory (“a.c.’’) component suffers reduction in 
magnitude and a change in phase. It is perhaps interesting to observe 
that for the large values of + the phase-shift approaches 7/2 asymptotically. 


§3. RELATIONSHIP wiTH NoisE-RepDuUCTION Factor, I™. 


Equation (2) is in itself of interest ; by comparison with (I), however, 


2 
where we found [2= ———, , we may show immediately 


1+ a?’ 
Ase AS 4/(1= 24), sin (aibeu), ca ee ee ren hee 


2, 
| ~=tan l=) : 
A/ \1—I? 


where A, is the signal after drift time, +. The variation of the signal 
components and of I? are plotted in fig. 1, particularly for comparison 
with the later work. 

Equation (3) may be interpreted as follows. Consider a wholly dis- 
ordered stream (full “ shot effect’, [2—1); then the noise current 3512 
may be considered as distributed in a frequency spectrum according to 
the well-known expression 


ST222el dios fk een ae 


After drifting for a time 7, according to (3), the magnitude will have 
diminished to 


Ey a2el (lon Jo") Gf. 5 suditve, = teak aera (ase) 


On the other hand, the d.c. component, I, which we have effectively 
considered as “ smooth” by distributing the fluctuation throughout the 
frequency spectrum in (4), will, by (1), have become partially disordered 
according to the law 


SPSL rah 7 ae ee aie 
The net sum of (4a) and (46) is thus simply 2eIdf. The irreversible 
maintenance of full shot effect (equation (4)) in the beam may thus be 


regarded as the “ detailed. balancing ” of these two mechanisms always 
proceeding at the same rate at any instant, When we consider the 
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problem, as in (I), of the beam with zero shot effect, then only one process 
can occur and we thus witness the so-called irreversible approach towards 
full shot effect *. We remark also that the signal-noise (power) ratio 
1—[? 

ee 

With the equivalence demonstrated above the calculation of I in 
certain problems is simplified. 


will evidently deteriorate following a law ~ 


§4. Lemma—Casn or Inrrran [7?240. 
Before continuing to a detailed analysis of the problems in view we 
state a simple lemma. Thus far we have assumed that J"?—0 initially ; 
we ask now how the noise will grow if initially 


SE =21T2df = (0<I2<1). 


Hips: 


Variation of signal components and /” for accelerated electron-stream. 


If I? , be the noise factor due to transit calculated for an initially smooth 
beam (T2=0), then the overall factor (say I”?) is given by 
Bee) eae a te a) aig ek Me ets <i apss ask cae D) 
A simple demonstration may be given by splitting the current, T say, 
into an initially ‘“‘ noisy’ component of magnitude IJ} and an initially 
“smooth ” one of magnitude I(1—I%) (ef. (1), p. 565). 


§5. [2 iw Varrous PoTentTiaAL DISTRIBUTIONS f. 
In (I) we assumed that the electron beam was injected into the drift 
region at relatively high energy, ®, with £2=0.. In practice the electron 
stream is derived from a structure (“ electron gun ’’) similar in principle 


* vy. d. Ziel’s communication (loc. cit.) is on somewhat similar lines to the 
foregoing. He suggests, however, that I” is derivable a prior from asa ner 
equivalent to (2) of this paper. We have reason to prefer to sss ee Hie 
analyses as of individual significance, demonstrating also, however ee. Lest 
assumptions (e. g. no space-charge interaction) the equivalence shown BBS a 

+ It should perhaps be remarked that we are here throughout oe in 
the dynamics of an electron stream—that is, in the convective part a ir + Wee: 
The experimental data used below for comparison refer indeed to such currents, 
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to a simple thermionic diode, of anode voltage ®/e (measured from the 
potential minimum, strictly speaking), having a hole in the anode for 
projection of the electron beam. The smoothing in density of the emitted 
current fluctuations occurs through the medium of the variable potential 
minimum which is present between cathode and anode when the current 
is limited by space-charge. At medium frequencies the overall transit 
time in the electron gun is so short that a quasi-steady state may be 
assumed in the ‘‘ classical’ analysis (cf. North (1940), Schottky and 
Spenke (1937)) ; consequently one need not consider what region of the 
“valve ’ is dynamically responsible for the smoothing process. At very 
high frequencies, however, the very existence of a high degree of smoothing 
in diodes and long beams (Kompfner et al. 1946; Kompfner 1947) 
suggests that the action must occur effectively within a very small fraction 
of the total diode transit time. Further, the general assumption that 
the ‘compensating ” stream of electrons released is of almost zero energy 
(relative to the space-charge minimum) suggests that essentially only a 
small region about the potential minimum is significant for the initial 
smoothing. 

_ As a first analytical model, then, on the basis of these considerations, 
we shall assume that the space-charge smoothing is effectively complete 
immediately beyond the potential minimum. We would emphasize 
that no theory as yet exists for the detailed space-charge interaction of 
electrons in a valve at very high frequencies ; we hope in the future to be 
able to make some progress in this direction. If a satisfactory detailed 
theory is forthcoming we should then be able to specify with greater 
certainty the “ initial conditions ’’ of our present problem. 

We turn thus to consider the current fluctuations in the electron gun 
itself. We measure potential, V, transit time, 7, and distance, 2, from 
the potential minimum. In line with the foregoing discussion, we assume 
first a smooth current “ injected’”’ at x=0; this corresponds to the 
assumption that the smoothing mechanism can be idealized, in the first 
place, to occur at x=0. 


(a): Linear potential variation *. 
Let V=Cx ; we first require the law of transit-time, ty, from V=0 at 


«—=0 to potential Vy at v=x, say, with variation of initial (thermal) 
energy, E. 


m\ da : 
Thus T= Je) ic ALY a. Wie) San ae ee (6) 
=7)(1—(E/eV,)"”), ae uate KE, 


where 7, is transit time for E—0, Equation (6) then replaces the equation 
for 7 appropriate in (I) and in §2 of this paper, namely . 
Ty=7)(1—EH/2®). 
* It is of course understood that a linear law could not exist fundamentally 


in the immediate neighbourhood of the minimum ; it is chosen simply as an 
obvious overall first approximation, ig 
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To calculate I? we rely on the equivalence demonstrated in $3 above. 
We find in this case that $,= A,/Ay is given by 


A ive) me ie r 00 ae < / 
5, sin wt [ eh cos (F) 22 dx-+-cos wt | e” sin (F) Bais, ax, 5 (8) 
“0 0 


where p= : (eV,/k6) ¥2, 


assuming that the electrons leave the minimum, as is well-known, with a 
thermal energy distribution given by equation (1) above. 
Thus 


* co foo) 


5. ==sin ot {1a | e* sin ax a +cos ol {a | 


«—** cos ax ax} (3) 
Jo 


0 


where ad=wr(k6/eV,)"? now. 


0/2 a = ; 
Ssin ot {1—ae-*"* | cat} 4-008 wt {*X eoes), Kaan Peet.) 


“0 


a/2 : 2 a? . 
Hence Pe=i—{(1—aer#/ é ir) + eet, es eae hha): 


The function | edt is tabulated (e. g. Jahnke and Emde (1933), p. 32), 


0 
or Lash Miller and Gordon (1931) tabulate «” | edt. One may also 


n 0 
usefully use expansion in powers of a, showing in particular that for 
small a 
a We 
2a? (1—- =, (10—3 Sees 


—0-215a2-L0-024a4+.... 
The relevant functions in (10) and (11a) are plotted in fig. 2, with a as 
abscissa. We observe now that a phase change of 7 in the signal occurs 
as T—> 00. 


In comparison with experiment, we may take the following values for a 
typical electron gun used at very high frequencies : 


fz (number of transit cycles) > 3 


V =4 x 10? volts 


la 


8-2 102 °K. .°. kG 8-6 x 10-2 volts. 
(2 


> 


4 . a / r oy, —2 NVQ ¢ 4 re. 
Hence, in this case, a=67 \/(2-15) x 10°? = 0-28, and therefore 


PASO Cra Wei ay sa merge) 
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Since experimental work based on the beam derived from such a gun 
(cf. Robinson and Kompfner 1950) indicates [* ~ 0-02 this theoretical 
result might be considered most satisfactory, suggesting that the con- 
ditions postulated (20 at a—0, etc.) could be regarded as providing 
an acceptable model for the fluctuational behaviour at high frequencies, 

Before accepting this indication, however, it appears essential to 
consider the effect of closer approximation to the true potential distribution 


> 


in the valve “ gun ”’. 


Fig. 2 
h 
ho 
84 
6 


Quadrature 


oe Component of Signal 


Variation of signal components and I? when V=Cz. 


(b) The Langmuir—Child Potential Law. 

Child (1911) and Langmuir (1913) showed that the current, I, in an 
idealized parallel-plane diode, limited solely by space-charge, initial 
cathode emission velocities being neglected, is given by the now classical 


expression 0 1 a var me 
On Tee ee ae : 
where d: anode-cathode distance. 
. V,,: anode-cathode potential difference. 


We adopt this diode as a second model for the potential-minimum/anode 
region in the electron gun structure, where the minimum is now considered 
as a virtual cathode at a=0. The law of potential variation with « 
corresponding to (13) is then 
Vesa, et ecce ae ree eae 
where C is a constant for the structure under a particular current. 
The transit-time law with initial energy E (cf. equation (6) above) 


is then ae ae 


where v=H/e, the equivalent initial voltage, 
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Simple substitution transforms (15) to 


2 fy\ V4 m eo 1(mi!? 4/ (c/v)) 
p= = > 
SEN 2607 ke 


which after expansion, integration and some reduction yields 


T,=T9{1—0-85(v/V,)¥440-167(v/V)—0-038(0/V,)?-+....}. . (16) 


sinh!” ¢ dd, 


An alternative, more perspicuous, approach is to split the range of 
integration in (15) at z=, where € is the value of x for V—v. The two 
integrals may then be expanded and we obtain 


Ta(L) & 7(0-28(v/Vo)* 4), 


Ty (2) SY To(1—1-13(v/V,)#4+0-167(v/V,)....). a Pet l 7) 


These expressions themselves seem not without interest in view particularly 
of a recent discussion by Pierce (1948) of transit-time variation 
with initial energy in a valve with a potential minimum due to space 
charge. He measures time, f, as we do, from the place of zero field 
(minimum of potential) and then uses the equation 


ay 4a Ga 


a (18) 


(his equation (43) 
to determine the transit time. 

Where I is the diode current, and ¢ the transit time of an electron to x. 
This equation, however, as an application of Poisson’s law—remem- 
bering 0V/dx= 0 at the minimum—is only valid where all electrons 
(producing I) have the same velocity and ¢ is the transit time of one of 
them and it is thus quite inappropriate in the present situation. If, indeed, 
I be supposed produced by a stream of electrons leaving the minimum 
with zero velocity and ¢ be the transit time of a“ perturbed ” electron with 
positive initial velocity then it is obvious that the error in (18) is gross 
indicating in fact a variation of dV/dx in the wrong sense. 

Equations (16) and (17) are also interesting in view of the common 
assumption (e.g. Rack (1938, p. 613)) that thermal energies produce 
negligible transit time variation in a valve except under extreme conditions. 
If, however, we set 

Vy ~ 10? volts ; v=k/e ~ 10-1 volts, 
then from (16): tp = 79(1—0-15), showing that in only 3 cycles transit time 
the difference would already amount to ~ 3 cycle. 

Returning to our analysis we proceed as in (a), using now 

; Ty AY T)(1—0°85(v/V_)"*), AERP (16a) 
in place of (6), 


remembering also (v/V))<1. In this case we find that the relative signal 


is now 

(va) oO 4 
S.=sin wt ( —a | ee sin ax ir) +-eos wt (« | e” cos ax da are) 
, 7“ Q 


“ 0 
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reai—4(1—aJ” «~** sin ax de) + («J « =" cos ax de) } (20) 
0 0 


where now a=wr 0°85(k6/eV 4) 4 
The integrals in (19) may aoe be expanded and we wae 


I'(7/4) 
3.91 Bs) 


and 


US pe 4 
S,=sin on(t = al i eres + cos ot (2G)4) .a— 
(21a) 


whence for small a 


~ 006602 = 0-002a5" Hey Maka ks | MESA rine ances. peermerte (2 Lea) 


S_ and I? could be calculated in fact as far as a=2 using such expansions 
but the work quickly becomes lengthy. For very large a, asymptotic 
expansions may be sought, but are not of great value here. It was, however, 
found on plotting accurately «~** that a very close approximation to this 
function can be found through the formula 


(jee 0<#<0°3 | 
=1—1-36(a—0-3)?, 0-3 <a <0-91 
=1-36(1-52—2)?, 0-91 <a <1-52 - + (22) 
=0. 1:52<a 


Thus an explicit integration of (19) and. (20) could be obtained : 


bre | | 5-44 
S,=sin wl< cos 0-91a (—0-011+- Ray (1—cos 0-61a) 


F 44 
+ cos ot sin ota ( 5 (l—cos 0-610)) | = apo) 


54.4. 2 
ies | { 0-011-4 e = (l—cos 0-61a)} 
a 


0: ea 29°57 


=1-+ (1—cos 0-61a)— 


(k= e0s.076 [2 om. pert one 


The values given by (23a, b) agree excellently with those derived from 
(21a, 6) in the common region calculated (a <2) and the relevant functions 
are plotted in fig. 3 

Turning to the numerical example quoted above, we now have a=2-175, 
whence from fig. 3, /’°—0-26. This figure is now greatly in excesss of the 
experimental value and it is further evident from the lemma of §4 that if 
I} 9 the situation is merely aggravated. 
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-(c) Discussion. P 
It appears clear that this disagreement arises at least in part from our 
assumption that the smoothed current “ emerges” from precisely 7=0 ; 
V=0. In fact it is evident that the dynamical smoothing process cannot 
be entirely localized, while a glance at equation (17) shows clearly the 
important contribution to the above value for I? played by the potential 
in the immediate neighbourhood of the minimum. This feature is 
evinced even more strongly when we recollect that with the true velocity 
spectrum the potential initially rises with a square law (cf. e. g. Langmuir 
and Compton 1931), since in that case if we again postulate a smooth 


current from x=0, the initial rise of !'* may be shown to be exceedingly 
rapid following a law 


I?~ wr log 


(24) 


yur? 
(y, = e° (where C is Euler’s constant), =1-781....). 


Fig. 3. 


ena . ea WT 8 
Variation of signal components and £ when V=Ca**, 


On the other hand in such a vital region it then becomes imperative to 
consider other factors in the physics of the fluctuation problem such as 
the Coulomb repulsion of the electrons ; at this point we recollect the 
lack of a detailed theory of space-charge interaction mentioned previously. 

Returning, however, to the Langmuir—Child law the EON: of 
this region may be readily assessed from fig. 4. The variation of ™ for 
the numerical values previously laid down has been computed when it is 
assumed that the smooth beam “ emerges ” from a potential value greater 
than zero up to “ thermal ” potential, h0/e. In particular we find that in 
the latter case, I” has fallen as low as ~ 0-006, well below the experimental 
value, restoring. some reasonable confidence in the postulat ed model. 
It is probably significant also to notice that the transition from bo 
quadratic potential law toa 4/3 power law occurs In the region (1+ 10) xk@/e. 
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It is clear that in the linear approximation of (a) small weight attached - 


to the very initial region and consequently the condition 1°00 at a=0 
was adequate. 
§6. CONCLUSION. 


The analysis presented above cannot be regarded as necessarily final in 
nature, but it is offered in conjunction with the models discussed as a 
contribution to the study of the behaviour of electron streams. 
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LXXX. Proportional Counters in a Magnetic Field: an Investigation 
of the Isomerism of Br® and a Measurement of the Fluorescence 
Yield of Krypton *. 
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Atomic Energy Research Establishment, Harwell +. 


[Received May 17, 1950]. 


SUMMARY. 


The proportional counter technique for measuring the energy of soft 
electrons has been extended to include higher energy radiations by placing 
the counter in a strong magnetic field. Suitable X-ray calibration 
sources consisting of long lived K-capture bodies which give lines at 
energies up to 50 keV. have been developed. 

In an application of the technique to a study of the isomerism of Br®, 
the conversion coefficient of the 37 keV. gamma ray was found to be 1-2, 
which is close to the calculated value for an electric dipole transition. 

In a second application, the fluorescence yield of the K shell of Krypton 
was measured. The value obtained was 0-67. 


§ 1. INTRODUCTION 


THE use of proportional counters for studying the energies of soft radiations 
has been described recently in a series of papers by Kirkwoood, Ponte- 
corvo and Hanna, and by Curran, Augus and Cockroft. For accurate 
measurements the great majority of ionizing particles must spend all 
their energy in the active volume of the counter ; if a particle escapes from 
the counter before coming to rest, the size of pulse it produces will not 
correspond to its energy and a distorted energy spectrum will be obtained. 
The limit to the energy which can be measured with a proportional counter 
is therefore set by the dimensions of the counter and the pressure of gas 
inside it. In this paper it will be shown that this limit may be extended to 
higher energies for any particular counter by placing it in a strong mag- 
netic field with its axis parallel to the field. The escape of particles through 
the side wall of the counter is reduced, since the track is coiled about an 
axis parallel to the field direction. By using a strong magnetic field, the 
diameter of a coiled track (for an electron moving perpendicular to the 
field) can be made much smaller than the range, as can be seen from Table I. : 
in which electron scattering is neglected. 


* Some of the results given here have already been reported briefly (Rothwell 


and West 1950.) 
+ Communicated by the Authors. 
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Of course, particles having a small velocity component perpendicular to 
the axis will escape through the ends of the counter even when the mag- 
netic field is strong enough to prevent escape from the side walls. The 
escape from the ends, however, is always much less important than escape 
through the side walls. The effect of Coulomb scattering along the elec- 
tron track complicates the situation. As far as escape through the side 
walls is concerned, the scattering will tend to make the lateral spread of 
any track somewhat greater than the radius of curvature given in the table 
and so will reduce the efficacy of the field slightly. Scattering will, how- 
ever, assist in reducing escape of particles along the axis since it prevents 
any particles from moving nearly parallel to the field for long and as 
a result being uninfluenced by the field. Considerable reduction in the 
number of electrons escaping is therefore to be expected when a 
magnetic field is applied. 


TaBLE I. 


Radius of curvature in 


Range in magnetic field of 

Electron air in 7500 gauss. 

energy NADP. (for an electron moving 

perpendicular to the 
field) 

50 keV. 4 cm. 0.1 cm. 

100 keV. 12 cm. 0-15 em. 
500 keV. 140 cm. 0-39 cm 


The possiblity that the multiplication properties of a counter might be 
affected by a strong magnetic field was considered. As discussed below it 
was found experimentally that the multiplication properties were only 
slightly affected by the magnetic field, in a way which did not invalidate 
the technique. 

The technique has been tested in conditions where, in absence of the 
field, the electron escape would prevent energy measurements, Electrons 
of energies up to 50 keV. have so far been studied with this technique in 
investigations of the nuclear isomerism of Br®° and of the Auger effect in 
krypton. 


§ 2. DESCRIPTION OF APPARATUS. 

The counters used in this work were made of soft glass with an aqua-dag 
cathode painted on the outside (Maze 1946). The wall thickness ae 
1 mm., the diameter 5 cm., and the active length was 34cm. The central 
wire was tungsten 0-005 in., in diameter. This type of counter was used 
in preference to the conventional metal walled counter to reduce the chance 
of active gas being adsorbed on the walls. Each counter was fitted with a 
glass tap and a ground glass joint by means of which it could readily be 
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connected to a filling apparatus for evacuation and refilling. Pulses were 
taken off the wire and fed through a head amplifier and a linear amplifier 
to a four channel electronic pulse analyser. (Cooke-Yarborough, Brad- 
well, Florida and Howells 1950.) The counter was placed in the hori- 
zontal gap of a large electromagnet. The pulses from the counter were 
fed through a long cable to the head amplifier which, with the rest of the 
electronic equipment, was removed as far as possible from the influence 
of stray magnetic fields. Several types of counter fillings were tried. To 
test their relative merits, a small quantity of A?’ was admitted to the 
counter at each filling and the pulse distribution from it was measured. 
The width of the peak in the distribution is very sensitive to the presence 
of impurities. The percentage half width at half height of the A? K 
capture peak was found to be 10:5 per cent by Hanna, Kirkwood and 
Pontecorvo. With argon and methane fillings, measured values of the 
half width were usually as high as 13-14 per cent, because the commercial 
methane used contained a considerable amount of impurities. Since 
methane cannot be purified very easily, carbon dioxide was tried instead. 
With argon and carbon dioxide fillings, values of the half width were 
reduced to 12-13 per cent; the carbon dioxide was purified by freezing, 
and pumping off any residual gas, repeating the process several times. 
When, in addition, the argon was purified by allowing it to circulate over 
heated calcium for several hours the half widths were reduced further in 
value to between 10-5 per cent, and 11-5 per cent. The best mixture was 
found to be about 10 parts of argon to 1 part of carbon dioxide ; too great 
a quantity of carbon dioxide gives a greater half width, too small a quantity 
causes instability of the multiplication factor. 


§3. Errect or MaGNetic FIELD ON THE COUNTER MULTIPLICATION 
PROPERTIES. 


It was thought that the influence of the magnetic field on the motion 
of the avalanche electrons and the electrons produced by the ionizing 
particle might affect the multiplication properties of the counter, so the 
following experiment was performed in order to test this point. A small 
quantity of A?’ was admitted into a counter, of 2 cm. diameter and 22 cm. 
active length, which was filled with a mixture of argon (21-8 cm. pressure) 
and methane (7-6 cm. pressure). The pulse distribution due to electrons 
emitted following the K capture process in A*’ has a peak at an energy of - 
2-8 keV. The range of these electrons is very small, so that virtually 
none escape from the counter. Any change in the pulse distribution when 
the counter is placed in a magnetic field must therefore indicate that the 
multiplication properties of the counter have been affected in some way. 
Fig. 1 shows the pulse size distributions measured with the counter 
axis parallel and perpendicular to a magnetic field of 7500 gauss and with- 
out the field. _ The points all lie on the same curve indicating that in this 
case the magnetic field does not influence the counter properties. Later 
in the work the experiment was repeated with a counter having a larger | 


o.P2 
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diameter, and in this case a slight reduction in pulse size was observed in 
the presence of the field. Fig. 2 shows the distributions obtained with 
this counter which had a diameter of 5 cm. and an active length of 34 cm. 
It was filled with a mixture of argon (24:4 cm. pressure) and carbon dioxide 
(2-4.em. pressure). The effect of the magnetic field is quite small and 
since the shape of the distribution is unaltered it is unimportant for many 
applications of the technique. 


Fig. 1. 
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Effect of a magnetic field (7,500 gauss) on the pulse size distribution from 
A®” in the 2cm. diameter counter. 


§4. Enercy CALIBRATION OF THE COUNTERS. 

Some means of calibrating the counters had to be devised before any 
measurements could be carried out. The usual method employs 
characteristic X-rays for this purpose. Since it is inconvenient to use 
an X-ray machine near a large magnet, long-lived radioactive isotopes 
which decay by K electron capture may be used as an alternative source 
of X-rays. A%’ (2-8 keV.) has been used by Pontecorvo, Kirkwood and 
Hanna, and He (5-9 keV.) by Maeder but X-ray lines of higher energies 
than these are needed. 
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K-capture isotopes which can be prepared by thermal neutron irradia- 
tion, relatively free from long period f activities, were used. Specific 
activities of a few milli-curies per gm. gave adequate intensity. Table IT. 


lists the K-capture activities which were obtained by irradiation in the 
Harwell pile. 


Fig. 2. 
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Effect of a magnetic field (7,000 gauss) on the pulse size distribution from 
A®’ in the 5cm. diameter counter. 


TaBce II. 
Half period Energy of Ka X-ray 
peouepe (days) (keV.) 
Cres 26-5 5-0 
Zn® 250 8-1 
Se 127 ae 
Pits 17 20-2 
Snt8 105 24-1 
Gdt8 155 41-6 
Yipte? 33 50:8 
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With the lower energy X-rays from Cr, Zn and Se, it was essential to use 
thin sources and to provide a thin window in the counter wall to admit the 
radiation. With thick sources, the X-rays emerge only from a thin surface 
layer and are masked by gamma rays accompanying the K-capture, 
process} which can emerge from the body of the source. Thin sources 
were prepared by dissolving the irradiated element in acid and then eva- 
porating a few drops of the solution on a mica foil. The sources were 
tested with a counter filled to a pressure of about 30 cm. of mercury with 
argon and carbon dioxide. 
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Pulse size distribution from Zn*®® source. 
Curve A: without critical absorber. 
Curve B: with critical absorber (-0005” Ni). 


Fig. 3 curve A shows the pulse size distribution obtained with a Zn® 
source (8-1 keV.). Curve B shows the distribution when a critical absorber 
(00005 in., nickel) which removes the Kf radiation is placed between the 
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source and the counter window. The Kf radiation has an energy of 8-9 
keV. and an initial intensity of about one eighth that of the K« radiation. 
As expected the Ka and K£ radiations are not resolved in curve A and the 
effect of the critical absorber is merely to reduce the width of the distri- 
bution slightly. The small peak on the low energy side of the main peak 
is due to the escape of the K radiation of argon (2-9 keV.) from the counter 
following the photoelectrce absorption of the incident X-rays in the K 
shell of argon. Similar distributions having well defined peaks are 
obtained from the Cr®! and Se™ sources. 


Fig. 4. 
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Pulse size distribution from Pd! source with, and without, a magnetic 
field (6,500 gauss). 


With the higher energy X-rays from Pd, Sn1¥8, Gd? and Yb!®, it is 
not necessary to use very thin sources and the radiations are sufficiently 
hard to penetrate the glass wall of the counter (1 mm. thickness) without 
suffering too much attenuation. 
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Background radiations, consisting of 8 rays or electrons emitted by the 
source, or electrons generated in the walls by quanta, give rise to a con- 
tinuous distribution of pulse sizes; these are greatly reduced by a mag- 
netic field. An electron which originates in the cylindrical wall is rapidly 
bent back into the wall by the field and spends only a very small amount of 
energy in the counter. Fig. 4 shows the distribution from a Pd1® 
source (20-2 keV.) with and without a magnetic field and demonstrates 


Fig. 5. 
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Pulse size distribution from Yb!% source with, and without, a magnetic 
field (4,100 gauss). 


this effect very clearly. A similar distribution is obtained from a Sn1% 
source. With the Gd!’ and Yb1® sources, the photoelectrons generated 
by the X-rays are of sufficiently high energy for a significant fraction of 
them to escape from the counter at the pressure used. Fig. 5 shows the 
pulse distribution obtained from a Yb1® source (50-8 keV.) with and with- 
out a magnetic field. It can be seen that the magnetic field greatly 
increases the intensity in the peak, as well as decreasing the background 
intensity. 
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§5. IsomERism oF Br®?, 

An investigation has been made of the well known case of isomerism in 
Br*?, by admitting a small quantity of carrier free methyl-bromide into a 
counter and measuring the pulse distribution from it. It is known that 
the 4-4 hr. isomer of Br8° decays into the lower state, an 18 min. f active 
body, with the emission of a 49 keV. y-ray and a 37 keV. y-ray in cascade. 
(Grinberg and Roussinow 1940, Berthelot 1944). The 49 keV. y-ray is. 
totally internally converted, but there is still some doubt about the value of 
the conversion coefficient of the 37 keV. y-ray and the type of transition. 
Berthelot has obtained a value of 0-64 for this conversion coefficient, by 
coincidence absorption measurements, and concludes that the transition is. 
magnetic dipole. Segre and Helmholtz (1949) suggest from parity 
considerations that it is more likely to be electric dipole, in which case the 
internal conversion coefficient should have a value different from that 
found by Berthelot. 

The method used for producing the active methyl bromide was suggested 
by Dr. Glueckauf (Glueckauf, Jacobi and Kitt, 1949). A solution of 
potassium bromide in acetic acid was irradiated for a short period in the 
Harwell pile. Argon was bubbled through the solution, carrying active: 
methyl bromide, generated during the irradiation, with it. The argon was. 
then passed through a solution of caustic potash and then over phos- 
phorus pentoxide to remove free bromine and water vapour. It then 
passed into a container cooled by liquid oxygen, where the methyl bromide 
was frozen down while the argon passed through. The container, still 
cooled by the liquid oxygen was transferred to the counter filling system 
and pumped free from air and other gases. The liquid oxygen was then 
removed and the methyl bromide allowed to evaporate into the counter. 

Curves A and B in fig. 6 show the pulse size distribution from the 
active methyl bromide in an argon filled counter (containing about 30 cm. 
pressure of argon and carbon dioxide), with and without a magnetic field 
of 7000 gauss. .This demonstrates very clearly the effect the magnetic 
field has in reducing electron escape. 

The five peaks in the distribution arise in the following way. 

‘the 49 keV. gamma ray is ‘totally internally converted, but the 
37 gamma ray only partially, there are peaks in the distribution at 
(37-+49)=86 keV., when both gamma rays are internally converted, and 
at 49keV., when the 37keV.gamma ray is unconverted and escapes from the 
counter. The chance of a 37 keV. y-ray being absorbed in the counter is 
very small and can be neglected in our counter for pressures of argon less 
than 1 atmosphere. Both the peaks are modified by the escape from the 
counter of a 12 keV. Bromine K X-ray which is emitted in some cases 
during the filling of the vacancy created in the K shell of Bromine when a 
gamma ray is internally converted.. Additional peaks will therefore be 
produced at (49— 12) =37 keV. andat (49—12+4-37)=(49+4-37 12)=74keV. 
when one Bromine K X-ray escapes, and at (49—12+37—12) when two 
Bromine K X-rays escape. Peaks are therefore expected at 37, 49, 62, 
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74 and 86 keV. The counter was calibrated with 20-2 keV. X-rays from 
a Pd? source, the position of the calibration peak is indicated in fig. 6 
and it can be seen that the five peaks occur at about the expected energies. 
Fig. 7 shows the pulse size distribution from the active methyl bromide in 
the same counter filled with krypton at a pressure of 66-5 cm. and carbon 
dioxide (4cm.). The escape probability of the bromine K X-ray is 99 per 
cent in the argon counter, but only 45 per cent in the krypton counter ; it 


Fig. 6. 
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Pulse size distribution from Br® in an argon filled counter. 
Curve A: with counter in a magnetic field of 7,000 gauss. 
Curve B: without magnetic field. 
Curve C: background from Br® (corrected for decay) with 
counter in a magnetic field of 7,000 gauss. 


can be seen that the relative intensities of all the ‘‘ K escapes ” peaks have 
been reduced in the distribution obtained with krypton. Conversion 
coefficients can be calculated from these curves. The 49 keV. y-ray is 
totally converted, which is verified by the absence of a peak at 
(37—12)=25 keV., so the conversion coefficient of the 37 keV. y-ray is 
given. by 

__ no. of converted 37 keV. y-rays _ Ig9+-I,4+1g, 

~ no.of unconverted 37 keV.y-rays  Iyg+I5, 


where I,. represents intensity in the 62 keV. (peak), ete. 


X39 


Proportional Counters in a Magnetic Field 883 


Table LIT. gives values of the conversion coefficient obtained in this way 
for various runs with several different counter fillings and bromine samples. 
In the case of the measurement with krypton a correction must be applied 
for reabsorption in the counter of the 37 keV. y-ray. The probability of 


escape is calculated to be 0-92 for a 37 keV. y-ray in the krypton filled 
counter. 


Tasre ITT. 
Counter filling Conversion coefficient 
A+CO, 1 3 
A+ OO, il “| 
A+CO, 1-2 
A+CO, 1-4 
Kr+Co, 1-2 
av. 1-2 
Fig. 7 
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Pulse size distribution from Br® in a krypton filled counter. 
Curve A: distribution in a magnetic field of 7000 gauss. 
Curve B: background distribution from Br® (corrected for decay) 
in a magnetic field of 7000 gauss. 


The main error in these determinations arises from the background, 
composed partly of 34 hour Br*, which can be measured when the 4-4 
hour activity has decayed, and partly of 18 min. Br® in equilibrium with 
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the 4-4 hour Br* activity, which cannot be separated from it. Curve C, 
fig. 6, represents the 34 hour background, corrected for decay. The con- 
version coefficients were evaluated 

(1) taking a curve close to C as the true background ; 

(2) taking a curve drawn close to the troughs between the peaks as the 

background. 

There is a difference of about 14 per cent between these two extreme 
values of the measured conversion coefficient, and it is the mean which has 
been tabulated in each case in Table ITI. 

The value of the conversion coefficient calculated for (a) a magnetic 
dipole transition, (6) an electric dipole transition, is 0-60 and 1-55 
respectively. (Hebb and Nelson 1940). Our value of 1-2 for the 
conversion coefficient proves that the transition is electric dipole rather 
than magnetic dipole. This conclusion is in contradiction with the 
results of Berthelot but agrees with a result of Lidofsky, Macklin and 
Wu (1950) which has just been briefly reported. 


§ 6. FLUORESCENCE YIELD OF KRYPTON. 


There are two different ways in which a vacancy in an inner shell of the 
electron structure of an atom may be filled; either a characteristic X-ray is 
emitted when an electron falls in from a higher level or the corresponding 
energy is shared between several electrons (Auger electrons), ejected from 
higher electronic levels. The fluorescence yield is defined for a vacancy in 
a particular shell of the atom in the following way. 


Fluorescent yield = Wx (y.....) 


__ No. of K, (L,M....) X-ray photons emitted 
a No. of vacancies in K, (L,M....) shell. ° 


Fluorescence yields were first measured by Auger in 1926. He investi- 
gated gaseous materials in a cloud chamber using X-rays to excite the 
atoms. When the incident X-ray is absorbed in the K shell a photo- 
electron of energy E,—Ex (where E,=energy of incident X-rays and 
EK, =binding energy of the K shell of the absorbing atom) is observed. It 
is accompanied in W, cases, 7. ¢. when a K X-ray is emitted, only by several 
low energy Auger electrons from vacancies in the L,M.... shells. (The 
fluorescence yield is much smaller for the L,M .... shells than for the K 
shell so that the probability of emission of one or more Auger electrons 
during the cascade process of filling a vacancy in the L,M .. .. shellsis very 
high). In the remaining (1—W,,) cases a moderately energetic electron of 
energy X,—E,y....(X,=energy of K radiation, Ey, .... —binding 
energyofL.M..... shells) is observed as well. The fluorescence yield was 
obtained from the ratio of the number of events of each type. 

A method for measuring fluorescence yields of solid materials was devised 
by Compton (1929). An ionization chamber is used to measure the 
number of X-rays (whose energy must be greater than the K binding 
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energy of the element being examined) absorbed in a block of the material ; 
the intensity of the fluorescent X-rays excited in the material is measured 
at the same time. After correcting for absorption of the incident X-rays in 
the L,M. ... shells the fluorescence yield is obtained directly. 

The proportional counter, as pointed out by Kirkwood, Pontecorvo and 
Hanna, provides an alternative method of measuring the fluorescence 
yields of gases. When an X-ray is absorbed by the K shell of the counter 
gas, and the vacancy caused by the emission of the photo-electron is filled 
with emission of Auger electrons, the full energy of the incident X-rays is 
spent in the counter. .When the vacancy is filled with emission of a 
K X-ray, however, there is a good chance that this X-ray will escape from 
the counter, provided the pressure in the counter is low. In this case an 
energy E,>—X, (where Ky is the energy of the incident X-rays and X, 
_ is the energy of the K radiation) is spent in the counter. Thus, two peaks 
are observed in the pulse size distribution and the fluorescence yield can 
be determined from their relative intensities. 

Fig. 8 shows a typical pulse size distribution obtained by irradiating a 
counter containing krypton and carbon dioxide (16 parts krypton: 1 part 
carbon dioxide) with characteristic X-rays (24:2 keV.) from a Sn!® K 
capture source. A magnetic field was applied to reduce the escape of 
photo-electrons from the counter and to reduce background radiations from 
the source. The energy of the K, X-ray of krypton is 12-7 keV. so when this 
radiation escapes a peak is obtained at 24:2—12-7=11-:5 keV. An 
unresolved peak due to the escape of the K, radiation of Kr (14-1 keV.) 
is also apparent. The K, radiation from the source was removed by a 
critical absorber (0-005 silver) placed between the source and the counter. 

Two corrections must be applied in order to evaluate the fluorescence 
yield, Wx, from the relative intensity of these two peaks. First, the 
intensity of the ‘‘escape”’ peak is reduced by reabsorption of the K X-rays 
of krypton in the counter gas. The probability of reabso ption was 
computed from the known absorption coefficients using unpublished 
calculations of B. Pontecorvo. This correction can be made small by using 
low pressures of gas and did not exceed 5 per cent in any of the measure- 
ments. A second correction was necessary to take account of the 
incident radiation absorbed in the L,M .... shells. Since the total energy of 
the incident X-ray is then spent in the counter, this process cannot 
be distinguished from K absorption without subsequent emission of a 
K X-ray. The ratio of the absorption in the K shell to that inthe L,M.... 
shells is known at the K absorption edge for any element and it is assumed 
that the ratio is the same for energies not much greater than the K binding 
energy. The ratio for Kr was taken to be 6-4. This value had to be calcu- 
lated from an empirical formula due to Jénsson (see Compton and Allison) 
since there are very few measurements of the X-ray absorption coefficient 
of krypton reported in the literature. It follows, therefore, that absorp- 
tion in the L,M.... shells constitutes a fraction 1/7-4 of the total absorption 
in krypton. A similar correction had to be applied in the technique 
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used with solid materials mentioned above. The possiblity of an error 
due to escape of energetic photo-electrons from the counter was checked 
by measurements taken with magnetic fields of 3000, 4600, and 5900 
gauss with a Sn!48 source. The relative intensities of the escape peak were 
0-55, 0-56 and 0-56 respectively, which shows that the escape is small. 


Fig. 8. 
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Pulse size distribution from Sn! source (24:2 keV) in a krypton filled counter 
(magnetic field 4300 gauss). 


Absorption of X-rays in the CO, was neglected since it is very small 
compared with that in krypton. 

The areas under the ‘‘ escape’ and “‘ normal ”’ peaks were measured 
by planimeter from graphs of the pulse size distribution. The total 
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number of counts recorded in each peak exceeded 6000 counts in every 
case. The results obtained with different pressures and using Pd13 (20-2 
keV.) and Sn118 (24-2 keV.) sources are given in Table IV. 

The mean value of the fluorescence yield from these results is 0°67. 
The probable error of the average is only 0-005, but the systematic error, 
which is mostly due to the L absorption correction, is certainly larger than 
this. 

Values given in the literature for the fluorescence yield of krypton are 
0-51 (Auger 1926) and 0-615 (Harms 1927). Locher (1932) gives a value 
of 0-585 obtained by interpolation from data for neighbouring elements. 
The value obtained by the proportional counter technique is therefore 
considerably higher than that obtained by previous measurements. 


TaB_e IV. 
Fraction of pulses in | 
Energy Deca gee “escape” peak aie 
é of HL k oe pulses | Corrected And fluore- 
peciuens ae ry Eee in for re- corrected | scence 
pare eae). | hem f0t es escape ”’| absorption | for L,M...... yield || 
[=eY’) mhonrury) peak of K {absorption 
radiation 
—hj 
| 
24-2 3,800 5-4 0-59 0-60 0-70 7 
24-2 3,000 7-4 0-55 0-56 0-65 | 
24-2 4,600 7:4 0-56 0:57 0-66 { 0-66 
24-2 5,900 7-4 0-56 0-57 0-66 | 
24-2 4,400 13-8 0:52 0-55 0-63 
20-2 3,800 2°9 0-59 0-60 0-69 | 
20-2 7,100 2-9 0-60 0-60 0-70 
20-2 3,800 5-4 0-57 0-59 0-68 0:68 || 
20-2 3,800 7-4 0-57 0-59 0-68 | 
20-2 4,400 13-8 0-55 0-58 0-67 


A measurement of the fluorescences yield of xenon was also attempted, 
but in this case the accuracy is greatly limited by background radiations. 
from the K capture source used. Fig. 9 shows the distribution from a 
Yb169 source (50-8 keV.) in a xenon filled counter. The main “ escape cn 
peak is at 50-8—29-7=21-1 keV. (since the energy of the K, ray of 
xenon is 29-7 keV.), and is much more intense than the full energy peak. 
The subsidiary peak at 16-4 keV. is due to escape of the K, radiation of 
xenon (34-4 keV.). This was shown to be resolved when the escape peak 
was investigated in greater detail than in fig. 9. ‘The K; radiation from 
the source, of energy 57-7 keV. cannot be distinguished in the main peak, 
but is easily seen in the “ escape ” peak, at i le 28-0 keV This. 
suggests that escape peaks may be useful in investigating radiations. 
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whose total energies are too close to be resolved directly. The remaining 
peak probably arises from some impurities in the Yb. The value of the 
fluorescence yield of xenon obtained was 0°81. In view of the background 
and the fact that only a single measurement was taken, the accuracy is 
much smaller than in the case of krypton. 


Fig. 9. 
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Pulse size distribution from Yb!® (50-8 keV) in a xenon filled counter (magnetic 
field 7,000 gauss). 
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LXXXI. The Eddy Viscosity in Turbulent Shear Flow. 
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SUMMARY. 


The comparative slowness of those processes in the turbulent wake 
which determine the turbulent energy balance makes it necessary to con- 
sider the whole wake flow in attempts to derive local values of such quan- 
tities as turbulent intensity and eddy viscosity. If the wake is in a state of 
moving equilibrium, 7. e. the development involves processes controlling the 
intensity and eddy viscosity to values appropriate to the particular stage 
of development, then the controlling mechanism must be of large scale both 
in space and time. It is suggested that the large slow eddies responsible for 
the phenomenon of intermittently turbulent flow form this mechanism, and 
a plausible control sequence is outlined in which the intensity of the large 
eddies controls the level of eddy viscosity. The detailed structure of these 
large eddies has been investigated by measuring the spectrum functions of 
the velocity fluctuation components in the wake of a circular cylinder. The 
group of large eddies is distinct, and it is possible to deduce that the large 
eddies are randomly placed in the mean stream direction, extend over most 
of one half of the wake, and have their vorticity roughly in the direction 
of maximum positive mean rate-of-strain. By a rough calculation it is 
shown that this type are near the state of energy equilibrium necessary for 
the operation of the suggested control process. Finally, the generation of 
Reynolds shear stresses by distortion of quasi-stationary turbulent fields 
is discussed, and possible reasons for the applicability of the eddy viscosity 
hypothesis are suggested. 


h 
§ 1. INTRODUCTION. 


Recent work by the author on the structure of the turbulent wake of a 
cylinder has shown that wake development is due to the existence of a group 
of large-scale eddies whose function is mainly convective (Townsend 1949 
and 1950). That is, they convect fully turbulent fluid from near the wake 
centre to the outer part of the wake, and simultaneously entrain previously 
undisturbed fluid bringing it into closer contact with the fully turbulent 
fluid of the wake core. The appearance of intermittently turbulent flow 
in the wake is a natural consequence of this process, and is the most striking 
manilestation of the presence of these eddies which contain only a 
relatively small fraction of the total turbulent energy. The complete 
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turbulent motion can be regarded as the superposition of a small-scale 
random motion containing most of the turbulent energy on this large-scale 
convective motion. Detailed measurements show that the intensity, 
scale and diffusive properties of the fully turbulent fluid are substantially 
constant at any section of the wake, irrespective of distance from the wake 
centre, and this is held to be a natural consequence of the mixing action of 
the large eddies, which causes the instantaneous properties of a sample of 
turbulent fluid to depend more on its past history of crossing and recrossing 
the wake than on its present position. This notion that the properties of 
the turbulent fluid are not directly controlled by the instantaneous shear is 
contrary to the usual assumptions of mixing-length theory, which are based 
on the notion of local and instantaneous equilibrium of the turbulent 
intensity. Ifthe statistical equilibrium of the wake properties is controlled 
more by the past history of the wake than by its present state, it is not 
possible to use this notion to determine the equilibrium turbulent intensity, 
and this paper describes an alternative approach, based on a consideration 
of the dynamic equilibrium of the large eddies. 


§ 2. NoTaTION FOR THE CYLINDER WAKE. 


The usual notation for turbulent flow with two-dimensional mean motion 
will be used, that is, the components of the mean velocity along the axes 
Oz, Oy, Oz are U, V, 0, and the instantaneous components of the fluid 
velocity are U+-uw, V+v, w. Fluid density is denoted by p, kinematic 
viscosity by v, and the instantaneous pressure by P-+-p, where p is the mean 
pressure. All mean values are taken with respect to time at a point fixed 
in space. 

For the cylinder wake, the coordinate system is such that Oz is along the 
axis of the cylinder, Oz is parallel to the direction of the mean stream, and 
Oy is perpendicular both to the direction of the mean stream and to the axis 
of the cylinder. Then, making the usual “ boundary layer” approxima- 
tions (Goldstein 1938), the exact equations of mean motion reduce to one 

CUma cure eee U 
—+—=y55- .-.-... . (21) 
ape IGG) Oy 
If it is assumed that the turbulent shear stress wv can be described in terms 
of a coefficient of eddy viscosity «, defined by 


a oU 

Uv —eE Oy 
then CL) = | 9 
gaye ae Mata ora (2) 


Using the same assumptions as in (1), the equation for the balance of 
turbulent energy is found to be 
Ww OR othe ere 0 te! ieee gay eb OP? 
WHI = (pee) t ws + 55, (ute tel) + OS =O, 
(3) 


3Q2 


892 A. A. Townsend on the 
where W is the mean viscous dissipation 
W=r(uV72u+oV20 + wv). 


If the flow is at a sufficiently high Reynolds number to permit the 
application of the theory of local similarity (Townsend 1948), then 


— Bane 
We 15(52) . Oe eer ee Ee 


and all the terms in equation (3) can be measured directly except (1/p) 
(dpv/dy), which can be obtained by difference. 


§ 3. Tae Mrxrnc-Lenetu Hypotuesis oF Locan ENerey EQuILIBRIUM. 


The equation for the turbulent energy balance (equation (3)) contains 
terms representing viscous dissipation, convection of energy by the mean 
flow, production of turbulent energy by working against the turbulent shear 
stress, energy transport by diffusive movements, and energy transport by 
flow down pressure gradients, and states that the energy density is constant 
in time as a result of a balance between the gains and losses of energy from. 
these processes. By assuming that the turbulence is so smail in scale that 
the equilibrium intensity is determined only by the local parameters of the 
flow, and further, that the motion is dynamically similar everywhere, it is. 
possible to derive the Prandtl relation between turbulent intensity and rate 
of shear, 
dU 
a 
This assumption of small scale implies that the terms representing 
convection of energy by the mean flow, energy transport by diffusive 
movements and energy transport by pressure flow are negligible compared. 
with those representing production of energy by shear and loss by viscous. 
dissipation. Ifthis assumption is not made, the energy equation produces no. 
definite results and the simplicity of the attack is lost. Several attempts. 
have been made to allow for the effects of convection and diffusion, but it 
has been shown experimentally that all the terms in the energy equation 
are of comparable magnitude, and that it is very unlikely that the mixing- 
length view of the energy balance is even a good approximation. This is. 
most clearly shown by considering the time necessary for a chance fluctu- 
ation in turbulent intensity to disappear. Since the local rate of production 
of _ turbulent energy is —wv(dU/dy) and the local energy density is. 
4 (u*+v*+w?*) this time must be of order 


a= Kit ee he Ne het oie nnn 


—} (wt e24-w*) / we, 
oy 


which experimentally is found to have a minimum value of 1-5 2/Us, i. e. 
more than the time necessary for the free stream to pass from the cylinder to, 
the point of observation (Townsend 1950). This estimate of the “ relax- 
ation time ” is supported by the experimental observation that the setting 
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up of dynamical similarity in wake flow is a very slow process, and is not 
complete even at a distance of 500 diameters downstream from the cylinder. 

With the knowledge that diffusive and pressure transport of energy are 
sufficiently large to change the energy balance by a considerable amount, 
it can only be concluded that energy equilibrium in the wake is a process 
involving the whole wake, and that the turbulent intensity in a sample of 
turbulent fluid is determined more by its past history of crossing and 
re-crossing the wake than by its present position. This conclusion, while it 
explains very naturally the observed invariance of the mean properties of 
the turbulent fluid across any wake section as a consequence of the statis- 
tically uniform history of all samples of fully turbulent fluid, makes the 
problem of computing the actual turbulent intensity from the energy 
equation almost impossible and a new approach is necessary. 


§ 4. THe Hypotuesis or Movine Equirprium In WakE Frow. 

While the kind of absolute local equilibrium considered in mixing-length 
theory does not exist in free turbulence, a moving equilibrium may be 
possible, 7. e. when the wake is in a steady state, the mean rate of production 
of turbulent energy over a whole section is controlled by the turbulent 
intensity in such a way that there is a stable equilibrium intensity for any 
wake section and that chance fluctuations from this intensity distribution 
will not persist as the wake develops. Considering the history of a small 
sample of turbulent fluid, it is possible that, when the turbulent intensity is 
above normal, increased loss of energy by viscous dissipation may outweigh 
the consequent rate of conversion of mean flow energy into turbulent 
energy and so tend to the equilibrium intensity sufficiently far down-stream. 
This may occur as the result of local changes in structure, similar to those 
implied in the Prandtl equation, or by overall changes caused by the 
increased turbulent intensity. The homogeneity of the fully turbulent 
fluid makes it very unlikely that local changes in structure control the equili- 
brium, but the group of large eddies could form a very effective large-scale 
controlling mechanism in spite of their comparatively small energy. These 
large eddies, conveying turbulent fluid from the wake centre into close 
contact with non-turbulent fluid from outside the wake, are primarily 
responsible for the lateral spreading of the wake and consequently control 
the rate of conversion of mean flow energy into turbulent energy. Taking 
the extreme examples, if they were to stop completely wake development 
would almost cease and the overall rate of conversion of mean flow energy 
to turbulent energy would become very small, while a large increase in their 
intensity would accelerate wake spread and increase the rate of turbulent 
energy production. This is easily shown by writing down the kinetic 
energy of the mean flow 


ioe) ile 
[°  4p(Uy—U)?dy=3pU8 
since 
( pU;,(Uy—U) dy=4pU2dC,, 
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where J is a parameter representing the wake width, d is the cylinder 
diameter, C’, is the drag coefficient of the cylinder. 

Then the total rate of production of turbulent energy over the whole wake 
section is 


d2 dl 
Us| kp(U)—U)2dy=4pU2 


072 dx’ 


proportional to the rate of increase of wake width. 


§ 5. Enercy BaLance oF THE LarGEe EDDIES. 


If the large eddies control the rate of production of turbulent energy, it 
follows that their equilibrium intensity is itself controlled, and it is simpler 
to consider the factors in this control process than the factors controlling the 
intensity of the fully turbulent fluid. It is clear that the large eddies must 
obtain energy directly from the mean flow, and it seems likely that the loss 
of energy by friction of such a large scale motion may be described by the 
same coefficient of eddy viscosity as is effective for the mean flow. Then the 
control process should lead to a balance in which the large eddies are main- 
tained in a condition of neutral equilbrium, the viscous dissipation only 
exceeding the energy gain from the mean flow by a sufficient amount to keep 
the wake in a condition of similarity. Any increase in the intensity of the 
large eddies leads to a too rapid spread of the wake, to increased production 
of turbulent energy, and finally to an increased eddy viscosity which 
decreases the intensity. If this is true, then for any particular mean 
velocity profile, the equilibrium condition will be of the form 


Uy!) 
€ 


aay 
=R, F 


where R,,’ is a pure number, characteristic of the shape of the mean velocity 


profile, u,, 1, are the velocity and length scales at the particular wake 
section. 


Applying the momentum condition for the cylinder wake, this becomes 


Us 


€ 


=R,. 


c 


This implies that the mean motion in a turbulent wake is always taking 
place at the effective Reynolds number R,,, whatever the cylinder Reynolds 
number, R,—U,d/v, may be, which is one way of regarding the observed 
invariance with cylinder Reynolds number of the properties of the mean 
flow. 

Experience with the laminar instability of a boundary layer suggests that 
the instability represented by the large eddies might be two-dimensional 
and periodic in space, but the experiments described below show that the 
large eddies are certainly three-dimensional and aperiodic. The 
implications of this superficially surprising results will be considered later. 
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§ 6. EXPERIMENTAL ARRANGEMENTS. 


The measurements were all made in the wake of a circular cylinder of 
diameter 0-159 em. placed in an air-stream of velocity 1280 em. sec. 
and usually at a position 500 cylinder diameters downstream. 
The air-stream was produced in the small wind-tunnel in the Cavendish 
Laboratory, of section 38cem.x38cm.. At the position of observation, 
the wake flow is very nearly in a state of similarity, and it is believed that 
the results obtained are substantially independent of the initial flow 
conditions close to the cylinder (Townsend 1950). 

Turbulent velocity fluctuations were observed with the hot-wire 
anemometer, in both the single-wire and X-wire forms. Since a straight 
hot-wire responds very nearly to the component of the velocity fluctuation 
at right angles to its length and in the plane of the hot-wire and the mean 
flow direction, a single hot-wire at right angles to the mean stream 
registers fluctuations of w, one in the direction (1/4/2, — 1/,/2, 0) 
registers fluctuations in (w+-v)/V 2, and one in the direction (1/4/2, 1 /+/2, 0) 
- fluctuations in (w—v)/V2. In this way the fluctuations components in the 
directions of the three axes, i.e. u, v, w, and also the components (u+v) [v2 
and (w—v)/V 2 could be detected and their intensities measured. 

To obtain the intensity spectrum of a fluctuating electric signal, several 
methods are available. For simplicity, a zero-beat heterodyne analyser was 
used, a circuit in which, by means of a balanced modulator, the amplified 
turbulence signal is multipled by a sinusoidal signal from a beat-frequency 
oscillator, and the intensity of the resultant measured after passage through 
a low-pass filter. Since only the sum and difference frequencies appear in 
the resultant, the effect is that the original turbulence spectrum is displaced 
along the frequency axis by the oscillator frequency, and the energy 
transmitted through the low-pass filter is proportional to the energy 
originally contained in a band of frequencies, of central frequency the 
oscillator frequency and of band-width twice the cut-off frequency of the 
filter. This form of spectrum analyser is believed to have advantages for 
the study of the low-frequency part of the spectrum. 


§ 7. PROPERTIES AND STRUCTURE OF THE LARGE Stow Eppigs. 


From previous work on the structure of the turbulent wake it is known 
that the large eddies, which are responsible for producing the intermittently 
turbulent flow in the cylinder wake, contain only about 5 per cent of the 
total turbulent energy and extend over the whole width of the wake. The 
spatial extent of the circulation in the down-stream direction appears to be 
even greater. Using circuits designed to measure intermittency factors 
directly, it has been shown that the occurrence of fully turbulent flowin one 
half of the wake is uncorrelated with its occurrence in the outer half (‘Town- 
send 1950), and it is likely that the eddies in the two halves are substantially 
independent. This result makes it unlikely that the large eddies represent 
a periodic instability similar to the Karman vortex street, but more definite 
evidence is found in spectrum measurements (figs, 1-4). 
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Spectrum functions of w-component. 


The spectrum of the uw-component of the velocity fluctuation has a 
characteristic form, clearly due to the superposition of two groups of turbu- 
lence of widely different scales (compare, for example, the spectrum of 
isotropic turbulence (fig. 5) which may be regarded as the spectrum of a 
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random collection of eddies with a continuous distribution of size), and it is 
natural to identity the low frequency component of low total energy with the 
group of large slow eddies previously postulated to explaim the experimental 
observation on diffusion and intermittency. Within the rather large 
experimental error, the spectrum function of the large slow eddies is a 
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Gaussian error function which is characteristic of an aggregate of similar 
eddies of simple structure, uncorrelated in position along the Ow axis, and it 
seems a fair inference that the large eddies are of this general character, 
7. e. are of simple structure and are distributed along the direction of flow 
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randomly and not in a periodic or quasi-periodic manner. The initial 
surmise that the large eddies might be similar to a vortex-street must be 
discarded, and it is necessary to study the eddies more closely. 


Fig. 5. 
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Further information can be obtained fromthe spectrum of the v-component 
of the velocity fluctuation. An assembly of randomly oriented simple 
eddies has a spectrum function for v which has a maximum at a non-zero 
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frequency, and this is confirmed both by the spectrum of the v-component 
in isotropic turbulence, and of the large wave-number components of wake 
spectra. However, the low wave-number components of the wake spectra are 
again similar inform to a Gaussian error function, and itis concluded that the 
large eddies must have a preferred orientation. This preferred orientation 
of the axis of vorticity cannot be in the yOz plane, or the spectrum maxima 
would be very pronounced, and, since the u-component of the large eddies is 
not small, the preferrred orientation is not along Ox. This suggests that 
one the of principal axes of rate-of-strain of the mean flow (which are at 
+45° to the Ox and Oy axes) may be the preferred direction, and this. 
hypothesis is confirmed by measuring the spectra of the velocity com- 
ponents in these directions and along Oz. Considerable intensity of the low 
wave-number component is found in the spectra of w, and (in the positive 


half of the wake) of the (u—v)/V2 velocity component, while the intensity is, 


much less for (w+-v)/V 2, 7. e. along the axis of positive rate-of-strain. The 
preferred direction for the vorticity of the large eddies is approximately the 
direction of maximum positive rate-of-strain, which suggests that the 
source of the energy of the large eddies is the stretching they undergo as 
they are convected by the mean stream. 

This examination of the wake spectra shows that the large eddies are of 
sumple structure, are randomly placed in the direction of the free stream, 
and have a preferred orientation roughly along the principal axis of positive 
rate-of-strain. By transforming the spectra into the longitudinal corre- 
lation function, it can also be shown that the average extent in the Ox 
direction of a large eddy is considerably larger than the width of the wake, 
and the typical large eddy must be rather elongated in this direction. 


$8. MmcHANISM OF THE SLow Eppigs. 

Accepting the structure of the large slow eddies that has been outlined 
above, it is easy to see that such eddies receive energy from the mean flow, 
and may use it partly in increasing their kinetic energy and partlyin working 
against the shear stresses caused by the smaller eddies of the main body 
of turbulence. This gain of energy cannot persist indefinitely for, after a 
sufficient time, the stretching action will cease and the large eddy will be 
dissipated (fig. 6). The typical history of a large eddy may be something 
like this. Initially, appreciable circulation occurs by chance in a circuit 
with normal inclined at a small angle to the Oy axis, and, as the circuit is 
convected by the mean flow, the energy increases as the area of the circuit 
decreases. In the early stages of development, it is likely that the dimen- 
sions of such an eddy are too great for the effect of eddy viscosity to be 
appreciable, but, as the circuit area decreases, the velocity gradients tend to 
increase and loss by eddy viscosity increases in importance. At the 
intermediate stage of development, it is probable that the energy gain from 
the mean flow and the energy loss by eddy viscosity are comparable in 
magnitude, and that the energy is nearly stationary. When the circuit 
area can no longer decrease, the eddy loses energy and disappears. 
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The viscosity control mechanism which has been described previously can 
easily be modified to fit into this scheme. When an eddy reaches an 
appreciable intensity, its convective action will lead to local spreading of the 
wake and to increased transfer of energy from the mean flow to the small- 
scaleturbulentmotion. Theresultingincreasein turbulent intensity and eddy 
viscosity reduces and finally stops the growth of the large eddy, and it may 
be supposed that, during the intermediate period of the eddy development, 
it is near a state of equilibrium such that the rate-of-change of total energy 
is small compared with the gain from the mean flow or the loss by eddy 
viscosity. The aggregate effect of a succession of these eddies is that the 
eddy viscosity is maintained continuously at a value such that the mean eddy 
is in this condition of dynamic equilibrium, the stabilizing process being : 
increase of energy of large eddy ——> increased rate of spread of wake ——> 
increased rate of production of small-scale turbulent energy ——> increasing 
eddy viscosity ——> increased energy loss from the large eddy. The time- 
scales for all these processes are comparable, and it seems very probable that 
a control process of this kind does operate in the wake. 


Fig. 6. 
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Development of large eddy. 


Confirmation of this notion can be obtained by inquiring whether eddies 
similar to the large eddies of the wake could be in a condition of dynamic 
equilibrium. Experimentally, the large eddies are found to be long 
cylindrical structures, generally similar to the eddy described by 


u=v=Az exp {—}o%(y2-+22)}, 

w= — Ay exp {—307(y?+2*) }, 
which is an eddy with vorticity roughly in the direction (1//2,— 1 [V2, 0) 
and unlimited in the Ox direction. In the absence of better information 
about the exact structure of a typical large eddy, let us consider the energy 
balance of such an eddy ina mean velocity field of constant shear dU/ay. 


Assuming the eddy viscosity « to be isotropic, the loss of energy by 
eddy viscosity is 


Da [| W dy dz, 
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nee) on (2) 8 
G0) @). 


and then \ 
D=3rpeA2a?, 


The rate of production of energy is 


= [-. pu dy dz 


If the energy is stationary, S=D, and the equilibrium equation is 


eU 

oe = Gaye 

aa Gea”. 
However, it may be argued that in a developing wake the equilibrium 
condition should be that the ratio of the eddy energy to the wake 
energy remains constant in the period of dynamic equilibrium. On this 
hypothesis, the equilibrium equation becomes 


OUe (ales ee 
ay tala |= 


where U,/x is the time taken for the mean stream to travel from the 
virtual origin of the wake to the point of observation. From observations 
on the wake at a cylinder Reynolds number of 1430, 


Care ecbed 


so the two forms of the equilibrium equation are not very different. 
At the same Reynolds number, it has been found that the mean value 
of the eddy viscosity is approximately 1/80 U,d (Townsend 1950). 
Under these conditions the equilibrium size of a large eddy is given by 


a2ad—=40 or 50, 


depending on which of the two definitions of dynamic equilibrium is 
used. In fig. 7, the function exp. (—xy’), representing roughly the 
distribution of energy in the large eddy, and the mean velocity defect 
U,—U are plotted on the same scale, and the equilibrium large eddy is 
seen to fit very well into one half of the wake. 
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This very crude calculation shows that eddies of the dimensions neces- 
sary to explain the slow mixing processes would be approximately in a 
state of dynamic equilibrium, and that the suggested viscosity control 
process is possible. In free turbulence, the occurrence of intermittently 
turbulent flow makes the existence of the controlling system of eddies 
especially obvious, but there is quite a possibility that similar processes 
may be at work in closed flows. The only experimental evidence for a double 
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structure of turbulence in a closed flow is the observation of MacPhail 
(1944) that, in the fully turbulent flow between rotating cylinders, there is 
superimposed on the small-scale turbulent motion a quasi-periodic velocity 
variation of scale similar to the original laminar instability described by 
G. I. Taylor (1923). This velocity variation probably controls the level of 
eddy viscosity in a similar way to the large eddies of the wake flow. The 
contrast between the periodicity of the large eddies in cylinder flow and the 
aperiodic nature of the wake eddies is probably due to the constancy in 
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time of the cylinder flow which allows the development of a periodic 
instability. In wake flow, before a periodic instability has developed 
appreciably, the wake dimensions have increased by too much to permit the 
continued growth of the original instability. 


§ 9. THE Concert or Eppy Viscosity. 


This description of the turbulent motion in the cylinder wake in terms 
of an eddy viscosity has the advantage that, assuming the dynamic 
equilibrium of the large eddies, the exact mechanism of the production of 
eddy viscosity is irrelevant, and, in principle, it is possible to estimate the 
equilibrium eddy viscosity without considering the energy balance of the 
turbulent motion. There is no strong experimental evidence against this 
use Of a coefficient of eddy viscosity to describe the small scale friction, for 
the scale and intensity of the small-scale turbulent motion are such that 
“time of flight’ and “mixing length” are small compared with times and 
lengths associated with the mean flow. What is certainly incorrect is 
the attempt to relate its value to the local parameters of the mean flow, 
implying local and instantaneous energy equilibrium, The evidence shows 
that the setting up of Reynolds stresses must be a fairly rapid process com- 
pared with the attainment of turbulent energy equilibrium, i.e. that 
transfer of energy from one velocity component to another under the action 
of the shearing must be rapid. The comparative slowness of energy transfer 
from the mean flow to the turbulence and then to heat suggests that 
viscosity is irrelevant to the process of stress production, and, for similar 
reasons, turbulent transfer between eddies of different sizes also must have 
little effect on the establishment of the shear stress. The problem then 
reduces to the distortion of a turbulent field. Ifthe axes of reference are 
parallel to the principal axes of strain, then the effect of an instantaneous 
distortion described by contradiction ratios J, m, n, is that 


E;=1f),  41=™N: C= ny » 


where lmn=1 and &€,, 7,, ¢, are the final vorticity components £ , 7, ¢o 
are the initial vorticity components. a 

For the wake, n=1, and if € is in the direction of positive strain, & 
increases while 7? decreases. In terms of the corresponding velocity 
components, ue decreases and v? increases, leaving w? unchanged, at least 
for small strains. This results in the production of a Reynolds shear 
stress, o(w2—v), set up instantaneously and proportional to the total strain, 
and thisis probably the first stagein the actual process. When the Reynolds 
stress developed by straining becomes large, the tendency of ordinary 
turbulence to isotropy should resist the tendency of the strain, to 
align the vorticity along the direction of position strain, and the 
eddy viscosity may represent the balance between the orientating action of 
the strain and the disorientating effect of the turbulent transfer between 
the velocity components. That the selective action of the strain acts on 
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the vorticity rather than on the velocity makes it clear why the turbulent 
motion in free turbulence is apparently “‘ isotropic.” The intensities of the 
velocity components along the conventional axes of reference are nearly 
equal, because these axes are at 45° to the principal axes of strain, and 
transfer between the vorticity components along the principal axes 
hardly affects the energy in the velocity components along the ordinary 
axes of reference. 

This notion that the eddy viscosity is a parameter describing a balance 
set up between the orientating effect of the shear and the disorientating 
effect of turbulent transfer is the last remaining possibility within the “ equi- 
librium” group of theories of the eddy viscosity (after energy equilbrium 
has been discarded), and its accuracy is dependent on a sufficiently rapid 
attainment of this equilibrium. Revelant experiments would be on the 
attainment of isotropy in a uniform field of, say, axisymmetric turbulence, 
but existing measurements are not sufficiently comprehensive and it is 
proposed to investigate this matter directly. If the tendency to isotropy 
is too sluggish to be an effective balance to the straining process, it may be 
necessary to discard the equilibrium concept of shear stress production 
altogether, and to suppose that, like the energy, the shear stress is deter- 
mined mostly by the previous history, probably by the total strain. While 
it would be premature to consider this hypothesis in detail, it should 
be mentioned that most of the conclusions concerning the mechanism of 
control of the large eddies would be equally valid, although the concept of 
eddy viscosity would then have less physical significance. 
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ABSTRACT. 


The effect of turbulence on ocean waves appears to call for an eddy 
viscosity which, while being large compared with the molecular viscosity, 
is small compared with the eddy viscosity found to apply to currents. A 
recent suggestion that the eddy viscosity should be taken as proportional 
to A> where A is the wavelength, seems to the author to be unsatisfactory. 
It is proposed, instead, that the coefficient of eddy viscosity N, applicable 
to waves, should be of the form N=Kca, where c is the speed of propa- 
gation, a is the amplitude and K is a constant. This form is indicated on 
dimensional grounds and is shown to be in conformity with v. Karman’s 
similarity hypothesis for shearing flow. With K of the order of 5x 10-5, 
the eddy viscosity would be large enough to account entirely for the 
observed rate of decay of ocean swell, while its effect on the attenuation of 
waves with depth would still be negligible. In the initial formation of 
waves, only the molecular viscosity need be considered. A possible effect 
of eddy viscosity on waves in the generating area would be to limit the 
steepness of the longer waves to a value less than the breaking steepness, 
even when the duration and fetch were unlimited. 


§ 1. [NTRoDUCTION, 

WuiteE molecular viscosity appears to be an important factor in the initial 
formation of surface waves on water, its dissipative effect on the long 
waves of ocean swell, for example, is negligibly small. If an eddy viscosity 
of the order of magnitude found to apply to ocean currents, on the other 
hand, were operative in the case of waves, their loss of energy would 
be much more rapid than is observed. Since it is well recognized in 
meteorology and oceanography that the eddy viscosity applicable to a 
particular motion depends on its scale, the suggestion has been made that 
waves may be affected by turbulence but that a smaller, and probably 
variable, eddy viscosity is applicable. 


§ 2, TURBULENCE AND Eppy VIscosiIty. 

The first quantitative expression of this idea would appear to be that 
made recently by Groen and Dorrestein (1950), who assumed an eddy 
viscosity proportional to \**, where A is the wavelength. This form was 
based on a result found by v. Weizsicker (1948) in the theory of locally 
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isotropic turbulence, and on an empirical law found by Richardson (1926) 
from diffusion experiments. V. Weizsicker was considering the transfer 
of energy from turbulence on a scale characterized by a linear dimension l 
to turbulence of scales smaller than J. By similarity considerations, he 
found the effective eddy viscosity to be proportional to 1*” for values of / 
within a certain range. Turbulence on a scale so small that dissipation 
by molecular viscosity is the dominating factor, or so large that transfer 
of energy from the mean motion to the turbulence predominates is excluded 
from the range to which the similarity law applies. In the case of the 
dissipation of waves, we are dealing with the loss of energy from the mean 
motion (7. e. the wave motion) and it is not clear, therefore, that the eddy 
viscosity law deduced by v. Weizsacker’s theory will apply to this problem. 

Richardson’s result, which was obtained originally for the atmosphere 
and has been found by Richardson and Stommel (1948) to apply also to 
the ocean in certain conditions, was that the coefficient of diffusion 
applicable to the dispersal of a group of particles is proportional to d*”, 
where d is the separation of a pair of particles. The process considered 
in this case can take place in turbulence which does not necessarily involve 
shear in the fluid and it is not apparent that the dissipation of energy in a 
problem involving shearing stresses can be deduced from it. 

A further difficulty is that, even if we assume an eddy viscosity propor- 
tional to I*°, the question arises : what is the characteristic length J for 
wave motion: Groen and Dorrestein have taken it to be A, but one might 
expect the amplitude a also to be a characteristic length for a given set of 
waves. 


§ 3. Wave MoTION AS AFFECTED BY MoLECULAR AND Eppy VIscostry. 


Let us consider two-dimensional surface waves in deep water. Let 
rectangular axes be taken, with Ow and Oy in the undisturbed sea surface 
and Oz vertically upwards. Let the waves be travelling in the direction 
of the negative x-axis, and U and W be the horizontal and vertical 
components of the particle velocity at any point. Then for irrotational 
waves of small amplitude 


U=— wae™ cos (wt+kzx), | 


W=-— ae sin (wt-+ka), (1) 


where w=2r7/T, T being the period, k=27/\, A being the wavelength and 
a is the amplitude of the surface elevation. ¢ is the time. 
The rate of dissipation of energy by molecular viscosity may be shown 
to be (Lamb 1932, p. 623) 
Ds 2h et ee ten) 7 ee eee) 


where y is the coetticient of viscosity and c=A/T is the velocity of 
propagation. If no external forces are acting, the amplitude of 
the waves will decay at a rate given by the factor exp (—2vkt), where 
v=p/p is the kinematic viscosity, p being the density. 
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To a first approximation, the components of the particle velocity are 
still given by (1), provided that 


k 
ely yg ys Ne nem) 


By inserting numerical values in (3), we find that @ is, in fact, negligibly 
small, except for very short waves. 

To consider the effect of eddy viscosity, we note that a wave motion 
is specified by the three quantities A, Tanda. We will assume, therefore, 
that if there exists an eddy viscosity N, applicable to wave motion, it 
is a function of A, T and a. If we assume a power law, N=K2a’TY 
where K is a non-dimensional constant, then, since the dimensions of 
N are L?T-1, we find 

a+t+p=2, y=—l1. 
Taking «=$=1 as the simplest solution including A and a, 


At 
N=K,7-=Kea. . 2. ef) 


It is suggested that equation (4) gives the appropriate form of eddy 
viscosity to apply to wave motion. We need not expect K to be an 
absolute constant but rather a non-dimensional parameter which varies 
only slowly compared with the other quantities involved. The above 
treatment is still useful if we can take K as constant over a limited range 
of the quantities A, a and T. 

The rate of dissipation of energy by eddy viscosity, replacing w by pN 
in equation (2), becomes, 

NaF nlosa ote. twee ae amr ear €3) 
The particle velocity components will still be given to a first approximation 
by (1), provided that the criterion (3) is satisfied with N replacing y, ¢. e. 


_Nk 


Sige eee te ae tb) 


It will be shown later that, with a value of K large enough for eddy 
viscosity to account entirely for the observed rate of decay of ocean 
swell, the condition (6) is still satisfied. 


§ 4. V. Karman’s Stminarity HyporHEsis APPLIED TO WAVE Morton. 


An alternative approach to the problem of the dissipation of wave 
motion by turbulence may be made by v. Karman’s similarity hypothesis 
for shearing flow, which has found successful application to turbulent 
flow through pipes and channels. Let w, v and w be the components of 
turbulent velocity and let square brackets [ ] denote values of functions 
of these components. Then if the mean velocity U is a function of z only, 
while V—W=0, the rate of dissipation of energy of the mean motion 
_by turbulence is given by 


dU 
= oO) —— : : * . ’ ° ° . al 
p=pluw] oe 2 a) 
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vy. Karman assumed dynamical similarity in the turbulence at various 
points in the field of motion and deduced that 


dU |dU 
= — | —— . . ° . ° . . 8 
ole dz\dz — (8) 
where : 
dU d2U 
2 9 
1=K,|=- / = (9) 


where K, is a non-dimensional constant. 
In the case of two-dimensional motion, in the xz-plane, (7) is replaced 
by (Lamb, loc. cit., p. 676) 


U oW oW 
p= {tury SE +) + to (G+ Se)pe 00) 


In applying this equation to wave motion, we assume that mean values 
of dU/dx etc., as well as [u?] etc. are taken over an element of volume of 
dimensions small compared with the wavelength and over a time small 
compared with the period of the wave motion. If U and W are given by 
equations (1), 


b= pkwae*{([u?]—[w?]) sin (wt+-kx)—2[uw] cos (wt+kx)}. . . (11) 


By analogy with (8), we should expect [uw] to be a function of dW/dx 


as well as of dU/dz. In this case, however, @W/dx—dU/d@z so that we may 
take 


,|dU| aU 
ea a Oz | Oz’ 
lou 2U| K 
where Lga=K, = / 12 ae 
where K, is a constant. 
Similarly, we will take 
dU | dU 
_ 1=0,+ Fl 55 [Fe 
K2| dW | dW 
ue Lev sales Eran cages 


where K, and K, are constants, introduced in a way similar to 
K,. C, and C, are positive constants, since [w?] and [w?] are essentially 
positive. From (11), we have, for the mean rate of loss of energy per 
volume from the wave motion at a given depth z, 


iW =f 
adW = =, dt 
T),” 
= pkw tate?" {(K3—K2) | sin (wt-+-ka) | sin? (t+ ke) 
-+2K7 | cos (wt-+-ka) | cos? (wt+-ka)} 
4 
ae pkeo*arette( Ke K2+2K%). 
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For the rate of dissipation per unit area of the sea surface, 
‘O : 4 
W=| _ dW de = — puta K3—K34+2K3). 2 . . (12) 
Hence | 
W=2pKk*c8a3 


as in equation (5) if 


04 
Ree (eR eRe. kk. 8) 


The equations of wave motion in the wz-plane, inthe presence of 
turbulence are (from Lamb, loc. cit., p. 678) 


dU Tada Deed we vee 
OE oye oe 
OW 1-0p 

ar ea 142 7 uwl=¢ 


where p is the wave pressure. To the ey approximation 

p= pgae™ cos (wt+ka) . 
The condition that equation (1) should still be valid, in the presence of 
turbulence, to a first approximation, is that terms such as d[uw]/dz 
should be small compared with terms ie as (1/p)dp/dx. 


We have = [uw]| = Kikw?a%e2** and sla gkaek*, Since K? is of 
the same order of magnitude as K, we may Se the condition as 
Kkw?ae?*? 
gkerrie 
and since kc?=g, this becomes 
Hager ele a ne. pn a ee) 
Apart from the factor exp kz, condition (14) is the same as that given 


by (6). 


We may therefore derive the equation for the dissipation of wave 
motion by turbulence from v. Karman’s similarity hypothesis without 
introducing the concept of eddy viscosity or the analogy with molecular 
viscosity. 

§ 5. THE Decay or OCEAN SWELL. 


5.1. Waves of a Single Period. 


We will consider waves of constant wave velocity c, travelling in the 
direction Oz and assume a quasi-stationary state, 7.¢. the wave energy 
varying with distance but constant in time at a given point. Then if E 
is the energy per area of the waves, the rate of loss of energy with distance, 
due to eddy viscosity, would be given by 
c dH 


sega 15 
5 ip (15) 
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assuming the energy to be propagated with the group velocity ¢/2. 
Since E=!pga? and taking W from equation (5), we have 
. da 4K ga? 

ic” mess 


using the relation kc?=g. The solution is 


i, a eee eee ee 


where a») is the amplitude at ~=0. In terms of the period T, since 
gF=2rnc, 
1 1 644K 
as a ce ee oe ae 
Ae ee (17) 
Sverdrup and Munk (1947), in their treatment of the growth and decay 
of waves, assumed that the effect of eddy viscosity was negligible and 
attributed the decay of swell to air resistance, the rate of dissipation per 
area being of the form 
Resse Oa, ts on Ge as tee 


where s is a coefficient and p’ the density of the air. If we replace W 
by R in equation (15), we are led to a solution ; 
A= OXPA Asp Gio jE} ia Mey 3) re 
or in terms of T 

a=, exp {—(4a-sp lpg l\x}... hs en 
A comparison of equations (17) and (20) shows the differences to be 
expected if the decay is due mainly to eddy viscosity or to air resistance. 
Both equations show that the decay is more rapid for waves of shorter 
period, the discrimination being greater in the case of eddy viscosity. 
In either case, therefore, if a spectrum of periods is present in the swell 
at a—0, the maximum amplitude occurring for a period Ty, then at a 
distance x, owing to the selective dissipation, the component of maximum 
amplitude will occur at a period T where T>T,. An increase of effective 
period with distance usually occurs with ocean swell, but Munk (1947) 
has shown that such an increase may be explained, in some cases at least, 

on kinematical grounds without any selective dissipation of energy. 


5.2. Hffective Period increasing with distance. 
Let us consider the general case of the period T, and thus the velocity c, 
increasing with x, without assigning the cause, and follow Sverdrup and 
Munk by assuming our dissipation equation to apply to the “ significant 


waves’. Ifc varies with x, equation (15) is replaced by 
1/ dE dc 
5(¢q, + BE) =—W. . . . . . . (21) 
Hence da lade 4Kg?a? 
=—5- (22) 
dix 2¢dz c4 


aT 
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Little is known of the way in which ¢ may be expected to increase with 2, 
but to obtain a solution we will assume that 

dc 

ene (a constant), so that c=c,+om. 


Then the solution of (22) may be shown to be 
1 1 SK g2x ( 1 1 ) 


: — — = = (| —— a 9 
Ge  aye,8? ~ T(é—Cy) \eg72 ce? 22) 


where dp, Cy are values at ~=0 and a,c are values at x. Putting gT=2re, 
1 __ 1, 1280tKe (1 1 
Ge Gees Tg Ty) ieee 172)” 


For dissipation by air resistance, it may be shown that the corresponding 
solution is 


(24) 


aT—aol'? exp {—(4rsp’/pgT)T)x}. -. . . . (25) 


5.3. Application to Observations of Swell. 


The above equations have been applied to the examples of swell given 
by Sverdrup and Munk (loc. cit., Table IV.). s and K have been computed 
from the recorded values of w, a), Ty, a and T, and are shown in Table I. 
For the examples of swell Nos. 1, 2, 3 and 5, equations (24) and (25) 
respectively were used. For Nos. 4 and 6, equations (17) and (20) were 
used, since T differs little from T,. The scatter of the values of s and K 
is about the same and one cannot decide on these data whether the decay 
law for air resistance or eddy viscosity fits the observations better. It 
is possible, of course, that both processes operate and are of the same 
order of magnitude. 

Taking the mean value of K, we would have 


INosote L0eCO Se) «et ne (20) 


As an example, for swell of wave velocity 10 m./sec. and amplitude 2 m., 
N=11-2 cm.2/sec., which does not appear unreasonable. Referring to the 
condition given by (6), we see that, with this value of K, the condition is 
easily satisfied. The influence of eddy viscosity on the rate of attenuation 
of the waves with depth is, therefore, negligible. 


$6. THE Tyrtra ForMATIoN or WavEs. 

Jeffreys (1925) showed that the early stages of the formation of waves 
could be explained if a transfer of energy from wind to waves took place 
through normal pressures on the wave profile. The rate of transfer of 
energy was taken to be of the form 


ESO Oa AUK ae) te We ol re ge 2 ea 
where V>c, s being the “ sheltering coefficient ’’ and V the wind speed. 


Jomparing A with D, the rate of dissipation by molecular viscosity, 
Jeffreys showed the existence of a critical value of V, below which waves 
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could not be formed. From observations he found V,,j,=110 cm./sec. 
(approximately), which corresponds to s=0-27. [If V>Vinin this treatment 
sets no limit to the height the waves may attain and they would grow, 
presumably, until they reached the breaking steepness of 20. 

From wind tunnel measurements on wave models of various ratios 
of 2a/A, Motzfeld (1937) deduced that A was of the form 


A=de'p'(ka}®2(V—c)%el. “yin Gb: eee 


TABLE I. 


Decrease in Height of Swell. 


2 a aw 5 6 7 8 
Decay | At end of fetch At distance x Air Eddy 
distance resist- visco- 
x Ampl. | Period | Ampl. | Period | ance s sity K 
Ay Ay a ail 
km. ‘ m. he m. sec. IOs lo 
2320 | 34 | 79 | O7 Tie clad 43 
1,480 2-4 8-1 FS 17-0 0-46 2-2 
1,950 3°45 8:7 1-0 17-0 1-39 Tel, 
1,110 2:35 10-5 1:35 Tee 1-20 6-1 
1,480 2°55 12-2 1-5 9-1 1-04 4-6 
185 1-25 9-0 1-05 9-5 1-84 9-3 
Mean 1-16 5-6 


The data in columns 1 to 6 are taken from Table IV. of Sverdrup and Munk’s 
paper (1947). Column 7 gives the value of s, computed from equation (20) 
or (25), assuming the decay to be due to air resistance only. Column 8 gives 
the value of K, computed from equation (17) or (24), assuming the decay to 
be due to eddy viscosity only. Ay 


with s’=0-014. Applying this result to the formation of waves, he found 
no critical wind speed in the region of 100 em./sec. but for any value of V, 
ka is a definite funetion of c, which is a maximum when c=V/3. For 
Ves 100 cm./sec., the maximum ka would be negligibly small, but would 
increase rapidly in the region just above 100 cm./sec. 


— 
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Neumann (1949), after considering observations of wind stress and 
wave steepness for various wind speeds, has proposed an equation for A 
of the form 


eee Dp WV C)r Ce se in, | he chen ef 2) 


with s’=0-095. Like Motzfeld, he found ka,,,,, to be a definite function 
of V and to increase rapidly when V was about 100 em./sec., but his 
values of ka were appreciably higher than Motzfeld’s. 

While differing in detail, therefore, all investigators agree that easily 
observable waves should first appear when V=100 cm./sec. approximately, 
and have a velocity c=V/3, i. e. 33 cm./sec. approximately. If we take 
a=0-1 cm., equation (26) would give N=1-85 x 10-4 cm.2/sec., 7. e. much 
less than the kinematic molecular viscosity v, which is 1-8 x 10-2 em.2/sec. 
We may conclude that eddy viscosity can play no part in the initial 
formation of waves and that only dissipation by molecular viscosity Dee. 
be considered. 

It seems unlikely that the relation N=Kca would hold down to values 
of N comparable with vy. We might, in fact, define a ‘“‘ Reynolds’ number 
for wave motion” by R,,=ca/v, and expect that turbulence would affect 
the wave motion only when R,,, was above a certain value. 


§ 7. THE GROWTH OF WAVES UNDER THE INFLUENCE OF WIND. 


Although Sverdrup and Munk have used an expression for A similar to 
Jeffreys’ (equation (27) for the transfer of energy from wind to waves 
by normal stresses for V>6 m./sec., they found it necessary to take a 
much lower value of s, i. e. s=0-013. They showed that transfer of energy 
may also take place by tangential stresses, according to the equation 


INGE ayes UTA et etal eo Maes meen 9) 


where y is the resistance coefficient and was taken as 2-6 x 10-°. Unlike 
A, A’ remains finite when c=V. Dissipation by molecular viscosity is 
negligible for waves for which c is several m./sec. If eddy viscosity is 
also negligible (as assumed by Sverdrup and Munk), the height attainable 
by waves under the influence of moderate and strong winds is limited only 
by the breaking steepness, provided the duration and fetch are both 
sufficiently great. In practice, it is observed that, while the shorter 
waves may be continually breaking, waves travelling with speeds 
approaching that of the wind seldom reach the breaking steepness, and 
usually do not exceed about 1/3 of that steepness. On Sverdrup and 
Munk’s theory, the explanation given is that sufficiently large durations 
and fetches do not, in fact, occur to enable the long waves to reach a 
oreater steepness. 

If eddy viscosity is the limiting factor, we see that the wave steepness 
cannot exceed an equilibrium value given by 


A+A’=W 
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Taking A from equation (27) and A’ from equation (30), we find 


PSB ye ae 
max 4p Kp? 


where B=c/V, and has been termed the “ wage age”. Thus if s and y? 
are constant, or are functions of 8 only, the equilibrium steepness 1s a 
function of the wave age, 8, only. This general result is still true if A is 
given by either of the alternative equations (28) and (29). 

From equation (31) the equilibrium steepness given by ka,,,x will be. 
less than the breaking steepness, for which ka=n/7, if 


PUL Cisse 
s(1—B)?+2y? ~ 4p 
inseiting the numerical values of p’ and p. We will consider two 
special cases of this condition : 
(a) If B=4 andy? is negligible, 
Ke 2:72 X10, 


With Jeffreys’ value, s=0-27, K>7:3x10-4. With Sverdrup and 
Munk’s value, s=0-013, K>3-5 «10-5, but this is hardly valid, since the 
value s=0-013 was derived empirically on the assumption that eddy 
viscosity was negligible. 

(Of 61, then 


ka 


ie, “33685010 


K>1-36 x 10-92, 


Taking y2=2-6 x 10-8, K>3-5 x 10-8. 

Comparing these inequalities with the value of K deduced in § 5.3, 
it seems possible that eddy viscosity may limit the steepness of waves 
travelling with a speed approaching that of the wind, but unlikely that 
it will limit the steepness of waves travelling at about one-third of that 
speed. 

A relation between wave steepness and wave age was, in fact, found 
empirically by Sverdrup and Munk. Until we have more data on the 
rate of transfer of energy from wind to waves, however, we cannot say 
whether the steepness of the longer waves is limited only by the duration 
and fetch, or by an eddy viscosity of the form we have been considering. 


§ 8 THE Errect on WAVES OF TURBULENCE DUE TO OTHER CAUSES. 


It has sometimes been assumed that the turbulence affecting waves 
should be attributed to the effects of wind driven or other currents. 
Alternatively, it has been suggested that the turbulence affecting currents 
is derived from wave motion. Neither assumption is made in the present 
treatment, the turbulence being regarded as inherent in the wave motion 
itself. The wave motion generates the turbulence and the turbulence 
reacts on the wave motion. The question still arises, however, whether, 
in the presence of a shearing current, the increased turbulence would not 
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greatly increase the effective eddy viscosity for waves. A clue to the 
answer may be found in equation (10), from which it is seen that for 
dissipation of energy by turbulence to occur, the mean values of terms 
such as [ww]dU/dz must be finite. This implies not only a correlation 
between w and w, but also a correlation between [ww] and dU/dz. It is 
assumed that, in the case of turbulence associated with the waves 
themselves, such correlations exist. 

If we now consider the turbulence associated with, for example, a 
wind-driven current, there will undoubtedly be some components in the 
spectrum of turbulence of scale comparable with, or smaller than, the 
wavelength of the waves. These components will contribute to [u*] 
and [w?] in equation (10) and may possibly have a finite value of [ww]. 
It is only, however, if the contributions to [w?]—[w?] and [ww] have 
components which vary in phase with dU/éx and 0U/dz respectively, that 
they will add anything to the dissipation of the waves. There is, therefore, 
no real discrepancy between the simultaneous existence of a large degree 
of turbulence due to other causes and a comparatively small eddy 
viscosity applicable to wave motion. 
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LXXXIII. Search for the Negative Proton in Fission, 


By J. M. C. Scorr* and E. W. TirrEeRTon, 
Atomic Energy, Research Establishment, Harwell, Berks f. 


[Received May 17, 1950] 


Durina the 1946 Cambridge Conference a paper was read by Broda, 
Feather and Wilkinson (1947) describing ionization chamber and chemical 
experiments designed to test the possibility of emission of a negative 
proton as an alternative to the emission of a B-particle and a neutron, which 
occurs in certain fission products and is responsible for the delayed 
neutrons. Recently Kuan-Han Sun (1949), who appears unaware of the 
experiments of Broda et al., has suggested that the delayed neutron 
emitter N17 might undergo negative proton decay. Still more recently, 
the situation has been reviewed by Feather (1949). The two competing 
reactions are : 


X$>YRii + (Baha) srs Qa) s yi eine) 
XPV Ls Po ee Qs ce, be 4G) on eth er oe 


If the negative proton has the same mass as the positive proton 
(Hi—m,c*), then from the known mass values (ef. Stephens 1949) 


Q.= Qe? 1:79 MeV atl ie en) 


Since Q, is positive for the delayed neutron emitted in fisson, it follows 
that reaction (2) is energetically possible and non-relativistic negative 
protons might be emitted. 

In the course of an investigation in this laboratory, involving the 
examination of large numbers of tracks formed by the slow neutron fission 
of U*®° in photographic plates, a search has been made for associated 
events which might furnish evidence of negative proton emission from 
a fragment. The delayed neutron emitters have lives ranging from about 
0-05 to 56 seconds (Hughes ef al., 1948) and, as the fragments come to rest in 
the photographic emulsion in a time less than 10~-'™ sec., it might be 
expected that the negative proton would be emitted after the nucleus has 
stopped. Ifa slow negative proton is emitted the question arises as to its 
ultimate fate. 

Among the possibilities are : 

(i) annihilation with a proton in ihe emulsion with emission of two 

y-Tays ; 
eee 
* Now at the Cavendish Laboratory, Cambridge. 
| Communicated by the Authors. 
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(ii) capture by a heavy nucleus such as Ag or Br and subsequent 
annihilation with a proton of the nucleus with a chance that one 
of the y-rays interacts with nucleons of the same nucleus leading 
to meson production and break-up of the nucleus to produce a 
pectats 

(ili) capture by a nucleus and conversion of the rest masses of the 
positive and negative protons into excitation energy of the nucleus 
in some other way, leading to a high energy nuclear explosion. 

(iv) radiative capture by a nucleus which, unlike the first three, does not 
involve the annihilation of the “ rest mass ”’ of the particle. 


The idea that negative protons may exist is derived ultimately from 
analogy with electrons. The hole theory, with its sea of occupied negative 
energy states, was first introduced as a complement to the Dirac electron 
theory and the discovery of the positron then fitted in naturally. The 
application of the same ideas to neutrinos is also generally accepted. In 
the case of protons, although the Dirac equation does not give the right 
magnetic moment, it nevertheless appears to be the proper way of 
describing particles with spin $. If the true relativistic theory of protons 
does include a hole theory and negative protons, it seems certain that they 
must be capable of annihilation with positive protons with emission of two 
y-rays, just like positrons. 

An objection to this view is as follows. If negative protons can be 
produced according to equation (2) which, omitting those nucleons not 
transformed, can be written : 


EONS Pee eet gh he ai(d) 


then it might appear that two neutrons in any nucleus can annihilate by 
way of virtual creation of a negative proton as an intermediate state 


NEENS PF Pema reat Mica 5 3 g5) 


and it follows that most known nuclei ought to explode. Thus if negative 
protons exist and are produced by Feather’s process (2), it would seem 
likely that they cannot annihilate with positive protons, and therefore that 
they do not bear the same relation to protons as positrons do to electrons. 
This argument would only rule out the creation by Feather’s process of a 
negative proton analogous to a positron. It is still possible to imagine a 
negative particle which cannot turn its rest energy into photons in this 
way. If created by Feather’s process, it would then be absorbed by 


process (iv), €. 7. 


brie e.0- energy i oo eee. (6) 
and gamma rays would be emitted similar to those from neutron capture, 
However, it does seem reasonable to postulate the existence of the anti- 
protons demanded by the hole theory, and meet the above objection by 
supposing that the neutrons occurring in nature can only react with anti- 
neutrons, and to assume that these (like positrons) do not normally occur 
in nature. 
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If we assume, therefore, that a negative proton can undergo annihilation 
with a proton, we may calculate the probability of the process by regarding 
it as the dropping of an ordinary proton into a hole. Even if the Dirac 
equation is not strictly valid for protons (in view of the magnetic moment 
difficulty) it seems to likely give the correct order of magnitude. The 
mathematical argument is the same as that for electrons, and gives a 
cross-section for negative protons (of velocity v<c) 


Ca Te IMAC, ah cla ek Beale CER 


and a lifetime of the order of 10-3 seconds in a material containing free 

protons. This time is much longer than that required for the particle to 
come to rest in the emulsion, so if the particle escapes from the fission 
fragment nucleus at all it ought not to be annihilated in flight. 

These arguments lead to the view that possibilities (i) and (iv) above 
would allow a negative proton to be observed in the emulsion as a track of 
appreciable length emanating from the end of a fission track. Possibilities 
(ii) and (iii) on the other hand would give rise to a track leaving the end of a 
fission fragment and terminating in a star of tracks, some of which would 
be at low energy. 

In the experiments which have been carried out in this Establishment 
some 6-5 x 10° fission tracks representing 1-3 x 10° fission fragments have 
been examined. About half were in [ford C, emulsions sensitive to 
protons of energy up to 40 MeV and half in D, emulsions which, under the 
conditions of the experiment, recorded protons up to a maximum energy 
of 2-5 MeV. No tracks longer than 3 (proton energy ~300 kV) were 
observed associated with the ends of fission tracks though many knocked-on 
protons (identified from the geometry of the events) were observed along 
the fragment tracks. One case of a 3-particle star having its origin at the 
end of a fission track was shown to be a chance juxtaposition of events. 
It was identified as being similar to many other stars found in the emulsions 
resulting from the decay of chains of «-emitters present as contaminant in 
the uranium solutions employed for impregnating the plates. 

From this result it appears that negative proton emission, if it occurs at 
all, is an improbable process. Accepting 0-8 per cent (Hughes et al., 1948) 
as the yield of delayed neutron emitters in the slow neutron fission of 
U*85 the experiment allows us to set an upper limit of 1 in 5x 10® to the 
probability ‘of negative proton emission relative to delayed neutron 
emission from these fragments. 


We thank the Director, A.H.R.E., for permission to publish this paper. 
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LXXXIV. The Emission of Short Range Alpha Particles from Light 
Elements under Proton Bombardment.—I11. The Reactions SLi(pa)®He 
and “Li(py)®Be*. 


By W. E. Burcuam and Joan M. FREEMAN, 
Cavendish Laboratory, Cambridge f. 


[Received June 14, 1950.] 


SUMMARY. 

The emission of short range alpha particles associated with the 
14:8 MeV. gamma radiation produced in the 7Li(py) reaction has been 
established. The energy distribution of these particles can be interpreted ~ 
as due to a level in the *Be nucleus at 3030-100 keV. energy, with a 
width of 2000 keV. 

Some observations on the energy release (3:97-0:03 MeV.), cross section 
and excitation function for the *Li(px)?He reaction are also presented. 


§1. INTRODUCTION. 
THE disintegration of “Li by protons is of particular interest since a well 
known resonance level of the compound nucleus *Be can easily be observed. 
The properties of this level have been studied in detail by Devons and 
Hine (1949) by observations on the capture reaction 
Wao Be* = SBe-p i ek ee (h) 
The energy of the gamma radiation emitted by the excited *Be* nucleus 
formed in this reaction was measured by Walker and McDaniel (1948) 
and components of 14-8 MeV. and 17:6 MeV. were found. The energy 
17-6 MeV. agrees with that expected if the *Be nucleus is left in its 
(unstable) ground state after the emission of the radiation. The 14:8 MeV. 
component could be due to a transition to an excited state of energy 
2-8 MeV. in the *Be nucleus, in accordance with the reaction 
7Li+1H > ®Be* > ®Be**#+y.. . . . « « + (2) 
Evidence from other reactions for the existence of such a state is sum- 
marized by Wheeler (1941), and a *Be nucleus left in this state is known 


to break up with the emission of two alpha particles of about 1-4 MeV. 
energy{. It follows that short range alpha particles should be associated 


+ Communicated by O. R. Frisch. 

+ (a) The energy of the alpha particle pair exceeds the energy of the excited 
state by 89-5 keV., which is the disintegration energy of the ground state of 
*Be according to the measurements of Tollestrup, Fowler and Lauritsen (1949). 

(b) Reaction (2) and the subsequent alpha particle emission might be referred 
to as the ’Li(py’,«)*He reaction, in accordance with the terminology adopted by 
Fowler, Lauritsen and Lauritsen (1948), but in this paper the abbreviation 
*Li(py)§Be* will be used. 


SER. 7, VOL. 41, NO. 320.—SEPT. 1950 38 


922 W. E. Burcham and Joan M. Freeman on the 


with the 14:8 MeV. gamma radiation observed in the Li(py) reaction, 
providing that this radiation results from a transition to the alpha-emitting- 
level of Be. In the present work, of which a preliminary account has. 
already been given (Burcham and Freeman 1949 a), such alpha particles. 
were found and their energy distribution was examined and compared 
with that predicted from the gamma ray spectrum. 

In the course of this work a number of check experiments were made 
on the disintegration of *Li by protons, mainly in order to assess the 
contamination by the lighter isotope of targets of separated ‘Li. The 
results of these observations are also presented. 


§2. EXPERIMENTAL METHOD. 
(a) General. 


The apparatus used in these experiments has been previously described’ 
(Burcham and Freeman 1949 b, 1950; subsequently referred to as I, IT). 
In some of the work the analysing magnet was turned on its side so that 
observations could be made in a horizontal plane and the target could be 
brought nearer the high voltage accelerating tube ; larger currents could. 
so be obtained. 


(b) Targets. 


Thin lithium targets were prepared by evaporating the metal in vacuo: 
on to copper discs. During transference from the evaporator to the 
target chamber these targets lost their metallic lustre and were presumably 
converted into lithium hydroxide or carbonate (cf. Bonner and Evans 
1948b). Targets of the 7Liand ®Li isotopes (similarly oxidized by exposure: 
to air) were also available. These were deposited on copper, silver or 
aluminium backings in the small electromagnetic separator at the Atomic 
Energy Research Establishment, Harwell, and had surface densities, 
calculated from the ion current, of between 10 and 100 ugm./cm.2 of 
lithium before exposure to the atmosphere. A fairly satisfactory thin 
target of “Li was also prepared by converting a 7Li sample into sulphate 
and evaporating 7 vacuo on to a copper button. 

Target thicknesses in keV. equivalent energy for protons and alpha- 
particles were calculated assuming a stopping power of 400 keV./mgm./em.2 
for lithium carbonate for protons of 400 keV. energy. This value is. 
consistent with the results of Warshaw (1949), with the figures given by 
Bonner and Evans (1948), and also with direct measurements made by 
comparison of the shape of the *Li(py) resonance curve with that obtained 
from a thick target of the same composition (cf. Tangen 1946). The 
targets used were of thickness equivalent to 10-90 keV. energy for protons. 


+ The surface densities given refer to the oxidized deposit except where it is 
otherwise stated. 
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__ §3. Resutts ror °Li. 
The disintegration of ®Li by protons leads to the emission of short range 
alpha particles and *He nuclei according to the reaction 


ene He ste | eS) 


Fig. 1. 


Particles per microcoulomb 


Particle energy x- mass/icharger’ (MeV) 


Alpha particles and *He nuclei from the reaction ®Li(pa«) He. Angle of observa- 
tion 135°. Proton energy (a) 330 keV., (6) 630 keV., (c) 930 keV. Target 
thickness 90ugm./cm.?. The shaded triangle gives the analyser resolution. 


Fig. 2. 


Particles per microcoulomb 


3 4 Ss 1-6 1-3 1-4 1S 6 17 


Particle energy X-mass/(charqey (MeV) 


Alpha particles and *He nuclei from the reaction ®Li(p«)*He. Angle of observa. 
tion 135°. Proton energy (a) 380 keV., (6) 450 keV. Target thickness 
15ugm./cem.?. The shaded triangle gives the analyser resolution. 


References to determinations of the energy release in this reaction and to. 
observations on the cross section, excitation function and angular 
distribution are given by Hornyak and Lauritsen (1948) and a later 
determination of the enegy release as 4-017-+--022 MeV. has been published 
by Tollestrup, Fowler and Lauritsen (1949). : 

Fig. 1 (a), (6), (c) shows the spectrum of particles obtained from an 
ordinary lithium target of surface density 90gm./cm.*. The curves 
were all taken with an analyser resolving power of 50 (4 per cent resolution 
in energy) and at an angle of observation of 135° with the proton beam. 
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Bombarding energies of 930, 630 and 330 keV. were used, and the *Li(pz) 
yield was monitored by simultaneous observation of the alpha-particles 
from the 7Li(pa)*He reaction. Fig. 2 (a), (b) gives curves taken at the 
same angle, and with the same resolution, for bombarding energies of 
380 and 450 keV. and with a target of surface density 15 pem./cm.? ; 
these two curves were used for detailed comparison with the results 
obtained with thin ‘Li targets. 

The two groups in figs. 1 and 2 are the *He and 3He particles emitted 
in reaction (1). The energy of these particles observed at 135° does not 
vary greatly for the range of bombarding energies 300-900 keV. This 
simplified the determination of an excitation function for the reaction 


Fig. 3. 


RRH 7He 90° 


(d) 


“He CG 


5°) 


FHe CG 


tHe 135° 


"Ls (pe) 


(e trom RRM 


Particles per microcoulomb (“Li reaction) 


300 400 500 600 700 800 900 Kefete) 


Proton energy (keV). 


Excitation function for the reaction ®Li(pa)*He ; « i 
pa)®He ; angle of observat ens 
(a) ?He nuclei, (>) *He nuclei, (c) 7He yield corrected for riba mate 
ae Ae *He yield corrected for motion of mass centre, (e) yield of 
i(~ax)4He reaction on a different scale ; © observati 
Roberts and Hafstad at 90°. iano oS 


over this range, since it was only necessary to set the analyser field for 
the peak of either group and to take counts for a series of bombardin 

energies ; checks were made at a number of points to ensure that ate 
analyser field was correctly set. The yield of reaction (1) so observed 
for the 90 ugm./cm.? target is shown in fig. 3; curves (a) (6) give the 
actual observations for the *He and ‘He particles, and curves (c) (d) are 
corresponding smoothed curves after correction for centre of mass motion 
The angle in the mass-centre system for which the observations ee 
actually made depends on the bombarding energy and varies from 139° 
to 143° over the energy range studied. The *He particles correspond to 
alpha particles emitted in the forward direction between ce of 30° 
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and 38°, and the difference in yields for proton energies above 500keV. 
suggests that the angular distribution of reaction (3) is asymmetric, 
though at lower energies there appears to be isotropy, in agreement with 
the results of Neuert (1938). The true asymmetry of the angular distri- 
bution is probably more than is suggested by fig. 3, since the curves of 
fig. 2, taken with a very thin target, show an even larger 4He/*He ratio at 
135°. This is due in part to the fact that the observed yield of a group of 
particles of given energy increases with target thickness as previously 
described (II, fig. 2), and the target thickness necessary for the maximum 
yield of the *He particles (of energy about 2-1 MeV.) is greater than that 
required for the maximum yield of the ‘He particles (energy 1-5 MeV.). 
In the case of the ®Li(px) reaction the target thicknesses necessary for 
maximum yield of the *He and ‘He particles (allowing for the difference in 
stopping power of the target for the two particles) are about 100 and 
60 wgm./em.? respectively for an analyser resolution of 4 per cent. If a 
rough correction for this effect is made for the 90 pgm./em.2 target a 
‘He/*He ratio in agreement with that obtained from the 15 pgm./em. 
target is found. 

The observations of Rumbaugh, Roberts and Hafstad (1938) on the yield 
of this reaction are also plotted (after fitting at the 430 keV. point) in fig. 3 ; 
their figures for the ?3He group at 90° show a much more rapid rise than the 
present results, although the increase in intensity of the “Li(px) reaction 
used for monitoring (fig. 3 (e)) was in agreement with their results. The 
resolution shown in the curves of fig. 1 is sufficient to make the relative 
excitation function reliable to about 10 per cent. As an extra check an 
excitation function was plotted using as detector an electron multiplier in 
which a mica screen was inserted to prevent the *He particles being 
detected. The *He particles were thus recorded with uniform efficiency 
over the energy range studied ; the excitation curve obtained agreed well 
with that shown in fig. 3. 

The energy release in reaction (3) was calculated from the curves of 
fig. 1, using ‘‘extrapolated’’ energies in order to eliminate the target 
thickness correction (cf. IT, §2 a). A special analyser calibration for 
“extrapolated” energies was used to replace that of I, fig. 8, and in this 
calibration the Q values for the reference reactions °Be(pa)®Li and 
*Be(pd)®Be were changed from the values given in I to the latest values 
(2-121 MeV. and 0-558 MeV.) published by Tollestrup, Fowler and Lauritsen 
(1949). A carbon deposit formed on the targets during bombardment and 
the thickness of this film was found as described in II, §2 (a) to be 8-3 keV. 
for protons ; a correction of 32-+12 keV. to the alpha particle energy and 
27+12 keV. to the *He energy has therefore been made. The results for 
the 90 ugm./cm.2 target are given in Table I.; no estimation of carbon 
thickness was made for the 15ugm./cm.? target, and energy releases 
deduced from the curves of fig. 2 have therefore not been included. The 
error in each of the stated Q values arises chiefly from the correction for 
carbon thickness and the uncertainties in proton energy (+10 keV.) and 
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in calibration of the analyser for extrapolated energies (+15 keV.) 
and amounts to +40 keV. Assuming that the observations provide three 
independent values for Q the mean value for the energy release in the 
SLi(pa) reaction is 
Q=3-:97--:03 MeV. 

This lies between the older values listed by Hornyak and Lauritsen (1948) 
and the latest figure given by Tollestrup, Fowler and Lauritsen (1949). 

The cross section for reaction (3) was obtained by comparing the yield 
of the *Li(pa) reaction with that of the 1°B(pa) and *Be(p«) reactions. 
The *Li target used had a surface density, estimated from the ion current 
in the electromagnetic separator, of 5-4 .gm./em.? of the isotope before 
oxidation, the !°B target was the 38 »gm./em.” target used in IT, and the 


TABLE I. 


Energy release in the ®Li(pa)?He reaction. 


Proton eey) 
energy 3He energy 4He energy Brom Resi 
(keV) (keV.) (keV) 3He energy *He energy 
322 2122 — 3-98 — 
622 2152 1482 4-00 3-95 
922 2162 1482 3-96 3°95 


*Be target was deposited by vacuum evaporation and found to have a 
surface density of 16 wgm./cm.2 by direct weighing. The yields from these 
: - N(E) dE 
targets were assumed to be proportional to the integral | puss taken 
over the observed alpha-particle group, N(E) being the number cf alpha 
particles recorded when the setting of the analyser corresponded to an 
energy E. The small variation of yield due to the slowing down of the 
incident protons in the target was neglected The integrals were evaluated 
graphically and it was found that, assuming isotropic distributions in the 
centre of mass system, 


oLi (px)(330 keV 


.) 
EAE ee Sere Oe 
SPB (pe(Ss0Kom 
"Li (p2)(330 keV.) _ 9 ga 
o°Be(pa)(5380keV.) = 


The cross section for the 1°B(p«) reaction for a proton energy of 330 keV. 
was taken to be 4-5 x 10-27 cm.? (II), and that for the *Be(pa) reaction for 
a proton energy of 530 keV. to be 1:3%10-%5em.2 (Thomas, Rubin, 
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Fowler and Lauritsen 1949). The corresponding values for the ®Li(pa«) 
cross section for a proton energy of 330 keV. are 4:9 10-26cem.2 and 
85x10 cm.*. The difference between the two values is probably due 
to systematic errors in the cross sections assumed ; the latter value exceeds 
the cross section of 5 x 10-?® em.2 obtained from the results of Bowersox 
(1939) and both are very much greater than the value of 8-5 10-2? em. 
deduced from the curves given by Rumbaugh, Roberts and Hafstad 
(1938). 


§4. REsuuts ror ‘Li. 
A search for the emission of short range alpha particles in reaction (2) 
was first made by bombarding a thin target of ordinary lithium with 
protons. Fig. 4 (a), (b) shows the results obtained with a target of surface 


Particles p¢r microcoulorhb (a) aso keV. 
Particles per microcoulomb (b) 390 keV, 


Particle energy X mass/(charge) (MeV) 


Alpha particles and *He nuclei from ordinary lithium bombarded by protons of 
energy (a) 450 keV. ; (6) 390 keV. The arrows show the position in which 
short range alpha particles emitted in the “Li(py)*Be* reaction would be 
expected for disintegration energies of *Be of 2900 and 3100 keV. 


density 17 ugm./em.2 (equivalent thickness about 20 keV. at the angle 
used), an angle of observation of 60° and an analyser resolution of 8 per 
cent ; curve (a) was taken at a bombarding voltage of 450 keV., at which 
the (py) resonance was observed, and curve (b) at 390 keV. for which there 
should have been very little effect due to the resonance level. The 
corrections for target thickness and carbon deposit were made by shifting 
the energy scale until the *He and ‘He peaks fitted the values expected 
from the energy release in reaction (3). The position of alpha particles 
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expected from the disintegration of a SBe nucleus with an energy release 
of 2900 and 3100 keV. is also shown and there is some slight evidence for 
the appearance of extra particles in this region on curve (a). The intensity 
of the incompletely resolved *Li(pax) group was estimated by integrating: 
over the distribution and was compared with that of the small group 
which might be due to the 7Li(py)*Be* reaction. Using the ®Li(p«) cross 
section given in §3, the thick target yield of disintegrations leading to the 
emission of 14-8 MeV. gamma radiation was estimated to be 2-1 10~*° 
per incident proton for a target of LiOH. The predicted figure is 
1:75 10-8, derived from the thick target yield of 5-35 10~° quanta 
(14:8-+17-6 MeV. radiation) per proton (Fowler, Lauritsen and Lauritsen 
1948), and the ratio 0-33 of the intensity of the 14-8 MeV. radiation to that 
of the total radiation at resonance (Walker and McDaniel 1948). The 
difference between the observed and predicted yields is probably due to 
the bad definition of the small alpha-particle group shown in fig. 4. 


Fig. 5. 
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Alpha particles from lithium isotopes bombarded with 440 keV. protons. Angle 
of observation 104°. (a) *Li; (6) Li. The shaded triangle gives the analyser 
resolution. Inset: Yield of gamma radiation and of alpha particles from 
*Li as a function of proton energy. 


In order to reduce the effect of the *Li(px) group, targets of the lithium 7 
isotope deposited on copper backings were bombarded. The lithium 7 
samples were enriched by a factor of about 500 over normal lithium, but 
owing to the low cross section for reaction (2) and the high cross section 
for the ®Li(px) process, and to the fact that the first targets used were more 
than 100 keV. thick, the yields from the two reactions were of the same 
order of magnitude. In the first experiments with thick targets, the 
background due to the *Li groups had to be subtracted, and the analyser 
resolution was therefore set at 12 per cent in order to increase the counting 
rate and improve the statistical accuracy. The curve given in fig. 5 (curve 
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(a)) which is that already published (Burcham and Freeman 1949 a, fig. 2), 
but transferred to an energy scale, is typical of these results ; a bombarding 
energy of 440 keV. and an angle of observation of 104° were used. There 
is clearly a group of particles centred at an energy of 1-38-40-08 MeV., 
and of width about 0:45 MeV., which cannot be due to contamination of. 
the target by *Li. The group observed from a thick *Li isotopic target 
is shown in fig. 5 curve (b) ; such contamination, in the correct proportion, 
could evidently account for the higher energy group shown in curve (a). 
In order further to establish the identity of the group of alpha particles 
of mean energy 1-38 MeV., the excitation function of this group was 
compared with that of the total gamma radiation over the range 350— 
700 keV. Part of the two curves is shown in fig. 5 (inset), and the close 
similarity between them definitely establishes the connection between the 
1-38 MeV. alpha particles and the resonant gamma radiation. With 
certain assumptions, discussed in §5, these results indicate the existence 
in *Be of a level, unstable with respect to alpha particle emission, at 
2900-200 keV. above the ground state and with a width of 900 keV. or 
less. The work described later in this section, in which the *Li contamin- 
ation was less important, modifies these conclusions. 

In order to see whether the higher energy group of fig. 5 (a) could also: 
be attributed in part to the “Li(py)*Be* reaction, the contamination of 
the *Li target by ®Li was estimated quantitatively, using the ®Li(d«) 
reaction as an indicator. The results obtained for the “Li target and a 
thin *Li target for comparison, were as follows : 


‘Li target: Intensity of *Li(dx) group=65+20 alpha particles/1000 nC.. 
‘Li target: Intensity of *Li(d«) group=203+10 alpha particles/1000 uC. 
SLi target : Intensity of *Li(px) group—115-+-5 alpha particles/1000 wC.. 


From these figures the intensity of the ®Li(p«) group on the ‘Li target was 
expected to be 37+11 per 1000 C., while the intensity found was 44+-5 ; 
the second group could thus, within the limits of error, be ascribed entirely 
to ®Li contamination. 

For a given isotopic enrichment the ®Li contamination becomes relatively 
less important in comparison with the effect from “Li as the target thickness 
decreases until a target of thickness equal to the width of the (py) resonance 
(14 keV.) is reached. Experiments were therefore made with thinner 
targets ; fig. 6 shows the results obtained with an evaporated ‘Li,SO, 
target. The thickness of the target is shown by the inset curve, which 
gives the gamma-ray yield as a function of bombarding voltage with a 
beam of 20 keV. spread in energy. Curves (a) (0) (c) (d) give the alpha- 
particle yields for bombarding voltages of 390, 450, 550 and 750 kV. 
taken with the analyser at an angle of 135° with the proton beam. The 
resonant character of the alpha-particle distribution and the presence of 
a group which may be ascribed to *Li contamination are again shown. 
The energy level in *Be deduced from these observations is centred at 
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3200 keV. above the ground state, and has a width of less than 1100 keV., 
though both these conclusions are uncertain owing to the presence of the 
6Li contamination. 


Fig. 6. 
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Alpha particles from thin *Li,SO, target bombarded by protons of energy (a) 
390 keV., (6) 450 keV., (c) 550 keV., (d) 750 kev. Angle of observation 
135°. Inset : Yield of gamma radiation as a function of proton energy. 
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Alpha particles from a thin target containing 50ugm/cm? of 7Li bombarded by 
protons of energy (a) 450 keV., (6) 390 keV. Inset: Yield of gamma 


radiation as a function of proton energy. The shaded triangle gives 
the analyser resolution. 
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At this stage more uniform depositis of 7Li of known density became 
available, and it was then possible to use thin targets with a higher 
percentage of lithium than in Li,SO,. Fig. 7 gives the alpha-particle 
spectrum obtained from an aluminium foil on which 50 ugm./em.2 of *Li 
had been deposited (as estimated from the ion current in the electro- 
magnetic separator). The inset curve gives the gamma ray yield, and the 
alpha-particle spectra were taken for bombarding voltages of 450 keV. 
(curve (a)) and 390 keV. (curve (b)) ; in both cases the angle of observation 
was 135° and the analyser resolution 4 per cent in energy. The yields for 
the 450 keV. observations were monitored by simultaneous observation 
of the total gamma-ray intensity. For a target of this thickness the 
curves of fig. 2 and the excitation function given in fig. 3 show that the 
yield of particles from ®Li contamination will change by only about 10 per 
cent between 390 and 450 keV., and the 390 keV. curve of fig. 7 has 
therefore been used as a background determination and subtracted from 
the 450 keV. alpha-particle spectrum. Since the resolution width of the 
analyser (60 keV. at 1500 keV.) is considerably less than the width of 
the alpha-particle distribution a normalizing factor for intensity must be 
used to give the true distribution. This factor was obtained by the method 
used by Thomas, Rubin, Fowler and Lauritsen (1949) in which the analyser 
was used to study particles scattered elastically at different energies from 
a thick copper target; comparison of the observed distribution with that 
expected from the Rutherford scattering law gives the normalizing factor. 

The ordinates of fig. 7, after correction for background, were each 
multiplied by the observed normalizing factor, and the resulting alpha- 
particle distribution is plotted in fig. 8 (curve (a)). The shape of this 
distribution may still be distorted at low energy, since the normalizing 
factor was obtained under conditions of high gain and low bias, while in 
the observation of the alpha-particle spectrum the bias was increased to » 
discriminate against scattered protons. Similar observations for proton 
energies of 450 keV. and 390 keV., and an angle of 135° were made with 
other thin targets containing 30-50 ugm./em.? of “Li; several of these 
targets showed no appreciable effect due to °Li contamination. Long 
counting periods were necessary with the analyser resolution used, and 
correction had to be made for the gradual growth of carbon deposit on. 
the targets during the experiments. This growth was estimated by periodic 
observation of the (yy) resonance curve. Fig. 8 (0) is the result of adding 
together all the thin target observations after correction for the carbon 
deposit, and after normalization and removal of background ; the curve 
has the same appearance as that of fig. 8 (a), although it is shifted in energy 
owing to the carbon correction. Most of the thin target observations 
showed the structure which appears in both fig. 8 (a) and fig. 8 (b), although 
the statistical errors (indicated in fig. 8 (b)) are large, and it would be 
permissible to draw a curve with a single peak through the points. The 
interpretation of this curve is discussed in §5; by integration over it the 
number of disintegrations leading to the production of 14-8 MeV. gamma 


\ 


932 W. E. Burcham and Joan M. Freeman on the 


radiation in the targets used (all of thickness greater than the resonance 
width) was found to be 1-05 10-® per proton, which is to be compared 
with the predicted figure of 1-75 10~® for a thick LiOH target. The 
difference is probably not significant in view of the uncertainty of target 
composition, and of the fact that the bias used in the alpha particle 
observations may have reduced the efficiency of detection below 100 per 
cent at low energies. The integrated yield was found to be about 20 times 
the peak yield with the analyser resolution used. 

A simultaneous observation of the alpha particle and gamma ray 
excitation functions was made for one of the 30 »gm./em.” targets for alpha 


Fig. 8. 


bel Tt a 
rat 


0.6 Ke) L4 1.8 22: 
ALPHA PARTICLE ENERGY (MEV.). 


ALPHA PARTICLE YIELD (NORMALISED) 


Alpha particles from thin 7Li targets bombarded by protons of energy 450 keV. 
Angle of observation 135°. (a) Observations shown in Fig. ‘7 (a), 
normalized and with background subtracted. (6) Combination of all 
normalized observations. 


particles of energy 970, 1235 and 1500 keV.; the curves, of which that 
for 1500. keV. is shown in fig. 9, were all very similar. 

The observation of resonant alpha particles at energies as high as 
1500 keV. makes it extremely unlikely that the effects observed can have 
been in any way due to protons multiply scattered within the analyser : 
the impluses from the 1500 keV. alpha particles were several times larger 
than those from the protons of 450 keV. 

During the course of this work the alpha particle spectrum was observed 
at a number of different angles, chiefly with the object, in the case of the 
thicker targets, of finding the angle at which the separation of the 7Li 
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and ®Li groups was best. Low resolution observations of this type showed 
no large differences in yield or spectrum of the resonant alpha particles 
over the range 45°-135°. 


§5. Discussion 


The results presented in §4 make it certain that there is a resonant 
emission of alpha particles associated with the 14-8 MeV. resonant gamma 
radiation from the “Li(py) process. Walker and McDaniel (1948) observed 
that the energy of the 14:8 MeV. radiation increased with increasing 
proton energy, and the emission of alpha particles must therefore follow 
that of the gamma radiation. The approximate agreement between the 
alpha particle and gamma ray yields makes it unlikely that there is any 
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7Li(py)*Be* : simulanteous excitation curves for emission of alpha particles and 
gamma radiation. Alpha particle energy 1500 keV. 


strongly competitive mode of break-up of the excited *Be nuclei left 
after the emission of the 14:8 MeV. radiation in accordance with 
reaction (2). Nghe 

The precise interpretation of the results given in fig. 8 (b) is difficult, 
since the alpha particles are emitted from a *Be nucleus whose momentum 
is partly determined by the recoil due to emission of a quantum of 14-8 
MeV. radiation. The angular distribution of the 14-8 MeV. radiation with 
respect to the incident proton is not known in detail, and the velocity 
distribution of the *Be nuclei produced in reaction (2) therefore cannot be 
calculated. On the assumption of an isotropic distribution of recoil 
momenta, the single alpha particle spectrum which would be observed 
at 135° when the ®Be nuclei are left with an excitation energy of 3000 keV. 
was calculated and found to be of width about 400 keV. The centre of 
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the distribution corresponded within 50 keV. with the energy expected if 
the gamma ray recoil had been neglected. In interpreting the alpha 
particle spectra obtained in this work, it has therefore been assumed, as a 
first approximation, that the gamma ray momentum may be neglected 
for the determination of the position of the energy level involved from the 
observed maximum yield. The deduction of the true width of the level 
or levels corresponding to the observed distribution of alpha particles is 
more difficult, since the relative importance of the recoil effect depends 
on the actual width of the level. If the width due to recoil (assumed 
constant over the spectrum) is 4 and the true width of the level is [) the 
observed width should be given approximately by 


P= V/(I?+ 4’). 


The validity of this formula was checked by graphical construction for a 
number of assumed values of I. 

The work of Walker and McDaniel indicates a level of *Be of width 
about 2100 keV. Neglecting the influence of recoil, in a first approxi- 
mation, this corresponds to a real width I, of 1050 keV. for the energy 
spectrum of single alpha particles. With 4400 keV. the observed 
width would be [=1120 keV., which is in good agreement with the value 
of 1100 keV. calculated from fig. 8 (b). The ®Be level reported by Walker 
and McDaniel is at an excitation of 2800 keV. above the ground state or. 
2900 keV. above the state representing two free alpha particles ; there is 
an uncertainty of about 200 keV. in this figure. Again neglecting recoil, 
the peak of the alpha-particle spectrum should be at 1450-+100 keV., 
which is in agreement with the value of 1560-+-50 keV. calculated from 
the mid-point of the curve of fig. 8 (b). The present results are therefore 
in agreement with the work of Walker and McDaniel within the limits of 
error of the experiments. 

The curve given in fig. 8 (b) shows evidence of a structure which was 
observed in a number of the thin target experiments made with an 
analyser resolution of 4 per cent. The statistical fluctuations were, 
however, always large, and it has not so far been possible to obtain reliable 
evidence for a splitting of the *Be level concerned. Further experiments, 
with larger intensities, and designed to display the total energy of the 
alpha particle pair, would be desirable to clarify this point +t. - The curve 


t (a) R. M. Littauer and P. Meyer have recently studied the alpha-particle 
spectrum in this laboratory using an analyser with increased resolution, solid 
angle, and coverage, and with large bombarding currents. Their results suggest 
a single level. 

(b) The possibility that the double peak shown in fig. 8 (b) is connected with a 
strong anisotropy in the emission of the 14:8 MeV gamma radiation is rendered 
unlikely by recent unpublished measurements by Devons and Lindsey on the 
angular distribution of this radiation. A study of the correlation between the 
direction of emission of the 14:8 MeV. gamma radiation and the following alpha. 
particles is being attempted. 
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of fig. 8 (b) is therefore taken to indicate a single level in °Be at an excita- 
tion energy of 3030-100 keV. above the ground state and a real width 
of about 2000 keV. 

A level in the nucleus *Be at an excitation energy of about 3000 keV.. 
above the ground state has been observed in the following reactions : 


11B(pa)®Be* Dee and Gilbert 1936. 

10B(da)®Be* Smith and Murrell 1939. 

4He(x«)*He Devons 1939. 

*Li(dn)&Be* Green and Gibson 1949. r 
2C(ya)*Be* Goward, Telegdi and Wilkins 1950. 


There is also evidence of a similar level in 8Be at nearly the same excitation.. 
which takes part in the break-up of *Li; this has been established by 
observations on the decay of *Li produced in the following reactions : 


*Li(dp)*Li Christy, Cohen, Fowler, Lauritsen and 
Lauritsen 1947 ; 
Bonner, Evans, Malich and Risser 1948 ; 
Baxter, Burcham and Paul 1950. 
MB(na)sLi Lawrance 1939. 
*Be(yp)sLi Titterton 1950. 


None of these reactions, with the possible exception of *Li(dn), suggests a 
level width of as much as 2000 keV.; most of them are consistent with a 
width of about 800 keV., and it is therefore not at present possible to 
identify the level concerned in the “Li(py)*Be* reaction with any previously 
known level of Be. It is, however, certain that this level lies within about 
200 keV. of those deduced from other experiments. In the particular 
case of the break-up of *Li, Baxter, Burcham and Paul (1950) using the 
same analyser as in the present experiments, found a level at 3130-50 keV.. 
which is to be compared with the value of 3030+100 keV. found in the 
present work. 

The excitation curves given in fig. 9 suggest that the alpha particle 
yield at resonance increases less sharply than that of the total gamma 
radiation. This might be due to an uncertainty in the background levels 
owing to the low counting rates, but it would be in agreement with the 
observation of Devons and Hine (1949) that the ratio of the 14-8 to the 
17-6 MeV. gamma radiation decreases at the 440 keV. proton resonance 
energy. 


ACKNOWLEDGMENTS 


We are grateful to Sir John Cockcroft and to Dr. W. D. Allen and 
members of his group at the Atomic Energy Research Establishment, 
Harwell, who provided the separated lithium isotopes. We also wish to 
thank Messrs W. Birtwhistle and D. D. Stewart for much technical 
assistance. 


936 On the Emission of Short Range Alpha Particles 
REFERENCES 


Baxter, A. S., Burcuam, W. E., and Paut, E. B., 1950, Phil. Mag. (to be 
published). j 

Bonner, T. W., and Evans, J. E., 1948, Phys. Rev., 73, 666. 

Bonnzr, T. W., Evans, J. E., Maticu, C. W., and Risszr, J. R., 1948, Phys. 
Rev., 73, 885. 

Bowersox, R. B. 1939, Phys., Rev., 55, 323. 

Burcuam, W. E., and Freeman, J. M., 1949 a, Phys. Rev., 75, 1756, 77, 287; 
1949 b, Phil. Mag., 40, 807 ; 1950, Ibid., 41, 337. 

Curisty, R. F., Conen, E. R., Fowzrer, W. A., Lauritsen, C. C., and 
LauRITsEN, T., 1947, Phys. Rev., 72, 698. 

Des, P. I., and GitpeErt, C. W., 1936, Proc. Roy. Soc., A, 154, 279. 

Devons, 8., 1939, Proc. Roy. Soc., A, 172, 564. 

Devons, S., and Hinz, M. G. N. 1949, Proc. Roy. Soc., A, 199, 56. 

Fowter, W. A., Lauritsen, C. C., and Lauritsen, T., 1948, Rev. Mod. Phys., 20, 
236. 

Gowarp, F. K., TeLeep1, V. and Witxyss, F. J., 1950. Proc. Phys. Soc., A, 63, 
402. 

GREEN, L. L., and Gisson, W. M., 1949, Proc. Phys. Soc., 62, 407 

Hornyak, W. F., and Lauritsen, T., 1948, Rev. Mod. Phys., 20, 191. 

LawrRance, A. M., 1939, Proc. Camb. Phil. Soc., 35, 304. 

Nevuert, H., 1938, Phys. Zeits., 38, 618. 

RumsaucGH, L. H., Roperts, R. B., and Harstap, L. R., 1938, Phys. Rev., 54, 
654. 

Smit, C. L., and MuRRE LL, E. B. M., 1939, Proc. Camb. Phil. Soc., 35, 298. 

TANGEN, R., 1946, Kgl. Norske Skrifter No. 1. 

Tuomas, R. G., Rusiy, 8., Fowier, W. A., and Lauritsen, C. C., 1949, Phys. 
Rev., 75, 1612. 

“Tirrerton, E. W., 1950, Nature, London., 165, 721. 

noah eis A. V., Fowier, W. A., and Lauritsen, C. C., 1949, Phys. Rev., 
76, 428. 

Waker, R. L., and McDanizt, B. D., 1948, Phys. Rev., 74, 315. 

WarsHAw, 8S. D., 1949, Phys. Rev., 76, 1759. 

‘WHEELER, J. A., 1941, Phys. Rev., 59, 27. 


LXXXYV. Short Range Alpha Particles from the Break-wp of ®Li. 


By A. 8. Baxter, W. E. Burcuam and E. B. Pavtt, 
Cavendish Laboratory, Cambridge t. 


[Received June 27, 1950.] | 


SUMMARY. 


~ 


_The energy spectrum of the alpha particles from the *LiBe* break-up 
has been studied using a magnetic analyser and thin targets. The most 
probable energy of the alpha particle pair emitted was found to be 
3150-L50 keV. 


§1. INTRODUCTION. 


THE nucleus °Li is an electron emitter of 0:89 sec. half-life which breake 
up in accordance with the reaction 


SLi > §Be*+6- > 24He+f- eer asa cg rea) 


with an energy release of 16-75 MeV. (Hornyak and Lauritsen 1950), 
The spectrum of alpha particles emitted in this process has been studied 
be several workers and an analysis of their results up to 1941 is given by 
Wheeler (1941 a). The main conclusion drawn by Wheeler is that the 
distribution of alpha particles at low energies can be accounted for by 
supposing the existence of a broad level in the *Be nucleus at an excitation 
energy of about 3 MeV.; reasons are given for believing this level to be 
distinct from another level in *Be at about the same excitation which is 
observed in the !B(px) reaction (Dee and Gilbert 1936) and in the 
scattering of alpha particles by helium (Devons 1939, Wheeler 1941 )). 
Recently new observations have been made on the energy distribution of 
both the alpha particles and beta particles emitted in reaction (1) ; 
Christy, Cohen, Fowler, Lauritsen and Lauritsen (1947) studied the 
momentum relations in the process and Bonner, Evans, Malich and Risser 
(1948) extended observations of the alpha particle spectrum to very low 
energies by projecting *Li nuclei into a cloud chamber. Cassels has 
obtained similar results, with better statistical accuracy, by projecting 
8Li nuclei into a photographic plate§. Hornyak and Lauritsen (1950) 
showed that existing data on the low energy alpha particle spectrum fitted 


+ Now at Chalk River Laboratories, Ontario, Canada. 
t Communicated by the Authors. 
§ We aregrateful to Dr. Cassels for informing us of his results before publication, 
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with their measurements of the beta particle energy distribution,-and that 
less than 2 per cent of all the B~ transitions from SLi proceed to the ground 
state of *Be. 

The present work was commenced early in 1948 with the object of 
improving the statistical accuracy of the existing observations of the low 
energy part of the alpha particle spectrum. The method of magnetic 
analysis used, in which the energy of only one member of the alpha particle 
pair from reaction (1) is observed, is inferior to methods which display the 
total energy release in each disintegration, and leads, as pointed out in 
§ 3, to a broadening of the spectrum. It is felt, however, that the present 
results are of some interest in showing the application of a magnetic analyser 
to a study of the *Li reaction. 


§2. EXPERIMENTAL. 

The alpha particle spectrum was investigated by means of the magnetic 
analyser described in other papers (Burcham and Freeman 1949, 
1950 a, 1950b). In the first experiments a lithium hydroxide target 
deposited on a thick copper backing was bombarded with deuterons of 
900 keV. energy. The energy release in the reaction ‘Li(dp)*Li is 
—0-193+0-008 MeV. (Strait and Buechner 1949, Paul 1950) and the ®Li 
nuclei are all projected forwards, the maximum angle of recoil being 34° 
with the deuteron beam. . The SLi nuclei, therefore, recoil into the target 
and target backing and a correction, amounting to about 100 keV. had 
to be made in consequence to the observed energy of the alpha particles 
emitted (cf. Bonner, Evans, Malich and Risser 1948). This correction is 
entirely eliminated in the cloud chamber and photographic plate experi- 
ments ; in later work with the magnetic analyser it was reduced by using 
thin foils of collodion as target backing. 

The collodion foil on which the lithium hydroxide was deposited was set 
at 34° with the deuteron beam and the radioactive alpha particles were 
observed in a direction at right’ angles to the foil. The thickness of 
both the foil and the target were measured by observing the 339 keV. 
19H(py) resonance with and without the foil in position, and were found 
to be 20-++5 keV. and 10+-5 keV. respectively for protons. The observa- 
tions were made using a suitable irradiation and counting cycle ; during 
the 2 second irradiation period the yield of the reaction was monitored by 
observation of the alpha particles from the ®Li(dx) reaction through a 
side window, and the radioactive alpha particles were counted, for different 
analyser settings, in a 2 second period following the interruption of the 
deuteron beam. It was found possible to bombard the collodion foils 
with currents of about 5 wA. without damage for a period of about an hour, 
and as the yield was large (about 400 radioactive alpha particles per second 
at the peak of the spectrum) good statistical accuracy was possible. The 
targets became blackened by deposition of carbon during the experiment, 
but repeat observations showed only a small shift of the spectrum over a 
period of an hour. The analyser resolution was 4 per cent in energy. 


EEE 
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~ §3. ReEsutrs. 


The alpha particle spectrum obtained is shown in fig. 1, in which the 
number of alpha particles for a given number of irradiation cycles, 
corrected by using the monitor counts to the same intensity of irradiation, 
is plotted against energy. The lower energy portion of the spectrum is 
shown on a larger scale. The total thickness of the backing foil and 
_ target was 30+-10 keV. for protons of 340 keV. energy, or approximately 
120+40 keV. for alpha particles of about 1500 keV., and in order to correct 


Fig. 1. 


Relative number of alpha particies 
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Alpha particle energy (MevV.). 


Alpha particles from the *Li-*Be* break-up; experimental points. The 
shaded triangle represents the analyser resolution. 


for this it was assumed that the average retardation of the alpha particles 
would be 60-+20 keV. Correction was also made for the varying trans- 
mission factor of the analyser over the spectrum by applying an experi- 
mentally determined normalizing factor (cf. Burcham and Freeman 1950 b). 
The effect of charge exchange on the low energy particles was also allowed 
for by using the.figures given by Briggs (1927) ; the ordinates at the very 
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low energies have been increased by the necessary factor and the corre- 
sponding high energy points reduced to allow for detection of singly 
charged particles. With these corrections and normalization the spectrum 
of fig. 2 (a) is obtained from the smoothed curve of fig. 1 ; this figure shows 
that the maximum yield occurs for an energy of 1575-425 keV. If the 
alpha particles are assumed to be emitted from *Be nuclei at rest the energy 
of the most probable alpha particle pair is 3150+50 keV. which is to be 
compared with the value of 3300 keV. found by Bonner, Evans, Malich 
and Risser (1948). The points obtained by these authors are plotted for 
comparison in fig. 2, with the assumption that the energy of the alpha 
particle pair divides equally between the two alpha particles, and it is 
clear that the width (about 1000 keV.) of the present distribution is 
considerably greater. This may be due to the corrections applied to the 


Fig. 2. 
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Alpha particle energy (Mev.). 


Alpha particles from the break-up of SLi; corrected distribution in energy. 
(a) present observations, (6) experimental points found by Bonner 
Evans, Malich and Risser (1948). The horizontal line shows the 
broadening expected in curve (@) as a result of target thickness and recoil. 


experimental points, although the increased width. must in part be due to 
the use of a target of thickness 120 keV. as a source of activity, and to the 
fact that the *Be* nucleus disintegrates into alpha paticles before it has 
lost its energy of recoil from the emission of the electron from 8Li. so that 
. e - : 
the two alpha particles in general have different energy. The recoil 
effect depends on the distribution of recoil momentum between the 8Be* 
nucleus and neutrino, and has been studied by Christy, Cohen, Fowler 
Lauritsen and Lauritsen (1947). With the assumption of the maximum 
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recoil energy of about 10 keV. and an isotropic distribution of recoiling 
*“Be* nuclei in a target of 120 keV. thickness, it was found by construction 
that a spectrum of the form given by Bonner, Evans, Malich and Risser 
(1948) (theoretical curve) would be broadened by about 200 keV. in E, 
without appreciable shift of the peak. The recoil effect is responsible for 
about 125 keV. of this broadening and the target thickness for about 
75 keV. The broadening is indicated in fig. 2, and correction for it 
brings the width of the present distribution into agreement with the cloud 
chamber work, although it makes it unsuitable for detailed comparison 
with a theoretical level formula. 

_ In deducing the position of the excited state of *Be* effective in the 
break-up of SLi, allowance must be made for the energy dependent factor 
(Q—2E,)° (Q=16-75 MeV.) which appears in the theoretical formula for 
the yield of alpha particles pairs of energy 2E, to take account of the 
varying probability of the 8-transition (Wheeler 1941). If the ordinates 
of the distribution given in fig. 2 are divided by this factor the peak is 
found to shift to 1610-425 keV., which corresponds to an energy for the 
excited state of 3130-50 keV., using the figure of 89+5 keV. for the 
energy of the ground state of *Be given by Tollestrup, Fowler and Lauritsen 
(1949). 
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LXXXVI. The Energy Release in the Reaction “Li(dp)*Li. 
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[Received June 27, 1950.] 


SUMMARY. 

In the reaction 7Li(dp)*Li the velocity of the centre of mass exceeds the 
velocity of the recoil nucleus in the centre of mass system and recoils will 
be observed only at angles less than a certain angle with the deuteron beam. 
This maximum angle has been observed for several deuteron energies 
and the Q-value for the reaction determined to be—0-187+0-010 MeV. 


$1. INTRODUCTION. 

THE energy release in the reaction *Li(dp)§Li has been determined by 
several methods. Rumbaugh, Roberts and Hafstad (1938) using a 
target enriched in the “Li isotope observed no protons of range greater - 
than 1:76 cm. with 0-86 MeV. deuterons. This sets an upper limit to 
the Q-value of 0-28 MeV. They also found no appreciable yield of SLi 
at deuteron energies less than 0-36 MeV. which sets a lower limit to the 
Q-value of —0-20 MeV. Strait and Buechner (1949) observed a proton 
group from bombardment of lithium with 1-36 MeV. deuterons. The 
protons were analysed magnetically and detected in photographic plates ; 
the Q-value reported for the “Li(dp)*Li reaction is —0-193--0-008 MeV. 

The present paper describes work done at the Cavendish Laboratory 
in 1948 employing a different method of finding the Q-value for this 
reaction, based on the observation of the maximum angle of recoil of the 
SLi nuclei. This method, which can be used to study a reaction for which 
the velocity of the centre of mass exceeds that of the recoil nucleus in the 
centre of mass system, was suggested by Newson (1935) and has been 
discussed by Bethe (1937) who shows that the maximum angle of recoil 
$y in the laboratory system is given by 


, M, 12 
S11 Op) =— ‘iat | (M+) ap +My |} Rk eee 
e i 


where My, M,, My, Mz are the masses of the initial nucleus and particle, 
and of the final particle and nucleus respectively and E, is the energy of 
the incident particle. For the 7Li(dp)8Li reaction the condition that ¢, 
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shall be real is that E,>|Q]| approximately, so that a maximum angle 
of recoil should be observable for deuteron energies exceeding 0-2 MeV. 
An examination of relation (1) shows that the sensitivity of the method is 
greatest when EH, is as small as possible and it may be shown that the 
accuracy of the Q-value so determined is least dependent on E, when Q is 
near zero. ‘The number of recoils per unit solid angle reaches a maximum 
at the angle ¢, for an isotropic distribution of recoils in the centre of mass 
system and this facilitates observation of 45. 


§2. EXPERIMENTAL. 


The experimental arrangement is shown in fig. 1. The deuteron beam 
passed through a grid which supported an aluminium foil approximately 
0-61 mgm./cm.* in thickness. A thin film of lithium oxide was deposited 
on the foil. SLi recoil nuclei produced in the target by deuteron bombard- 
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Apparatus used for detection of *Li recoils produced in the *Li(dp)8Li reaction. 


ment were received on a zinc sulphide screen set at an angle near the 
maximum expected. A movable shield in front of the fluorescent screen 
varied the angular range over which recoils could reach the screen from the 
target. The alpha particles from the break-up of the *Li nuclei received 
by the fluorescent screen produced scintillations which were detected by 
an RCA 931A electron multiplier placed close behind the screen. 
Periodic irradiations of 2-3 sec. were made and the alpha-particle 
activity on the screen was observed in periods of 3-2 sec. following each 
irradiation. During these counting periods the deuteron beam was 
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interrupted and a rotating shutter automatically moved between the 
target and the screen in order to prevent detection of any alpha particles 
directly from the target. 

The angle fixed by the edge of the shield was measured by replacing the 
target by a sheet of tracing paper illuminated from behind. The distance 
from the point at which the beam intersected the plane of the fluorescent 
screen to the edge of the shadow on the screen was measured with dividers. 

The movable shield allowed *Li recoils to be recorded at angles greater 
than that fixed by the edge of the shield but the exposed screen area was 
different for each setting of the shield. This correction was found by 
placing a polonium source in the target position and counting for different 


settings of the shield. The correction curve so obtained was used to allow 


_ for the change in total solid angle with angle of recoil. 

The air equivalent of the target material and supporting foil was 
measured before and after bombardment using a polonium source and an 
ionization chamber. The difference between the residual extrapolated range 
of the polonium alpha particles when the target foil was in place and 
removed was 0:81-+0-04 cm. air. 

A proportional counter viewed the .target continuously through a 
mica window in the side of the vacuum box. The counter was arranged 
to detect the alpha particles from the reaction *Li(d«)4He and served to 
monitor the conditions of irradiation. 

About 40 irradiations were performed for each setting of the movable 
shield. The background due to dark current in the multiplier and gamma 
radiation from the target was measured when the shield completely 
covered the screen. Runs were made with primary deuteron energies of 
950, 900 and 850 keV. ‘The energies of these deuteron beams after having 
-passed through the aluminium foil and lithium oxide target were taken to 
be 550+20 keV., 460-+30keV., and 400-+-30 keV. respectively. ®Li 
nuclei could escape from the target from a depth of about 0-4 mm. air 
equivalent, that is, about 10 per cent of the thickness of the lithium target. 
Thus, for example, deuterons of initial energy 950 keV. produced 8Li 
nuclei which emerged from the target only for a deuteron energy of between 
570 and 530 keV. This was assumed to give the uncertainty in E, since 
other sources of error were estimated to be less important. Scattering of 


the *Li nuclei before emerging from the target was neglected, although 


such scattering probably contributed to the uncertainty in the measure- 
ment of the maximum angle of recoil. 

The background count per irradiation was subtracted from the total 
count per irradiation and correction was also made for variations in the 
count from the monitor counter. Finally the count of delayed alpha 
particles for each angle setting was multiplied by the corresponding 


geometrical factor as determined with a polonium source in place of the — 


target. This gave the mean count per unit solid angle per irradiation for 
each setting of the shield. A check made of the period of the activity 
being observed gave a half-life of 1-3--0-5 sec. 


————— ee 
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§3. RESULTS. 

The corrected results of the various runs are presented in figs. 2 and 3. 
The ordinates correspond to counts per irradiation for 100 counts per 
irradiation in the monitor counter for a solid angle equal to that subtended 
by the screen in the fully exposed position. The values of the maximum 
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TABLE. 


Run | E, (keV.) | ¢p (degrees) Q (MeV.) 


460-+30 | 26-6-+1-5 | —0-193-40-024 
400-430 | 27-2+0-5 | —0-162-+40-013 
550-+20 | 27-840-8 | —0-215-+0-014 
460-430 | 27-7-40-8 | —0-182+0-015 
400-30 | 25-2+1-4 | —0-182-+0-021 


Ou 09 bo 


angle are listed in the Table together with the Q-value calculated for the 
reaction, ; ; 
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Fig. 3. 
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These results have been combined and the weighted average for the 
energy release in the reaction ‘Li(dp)*Li is found to be 


Q=—0-187+0-010 MeV. 


This value is in satisfactory agreement with the values of —0-20-+-0-03 
MeV. found by Rumbaugh, Roberts and Hafstad (1938) and of 
—0-193-+-0-008 found by Strait and Buechner (1949). 
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ABSTRACT. 


The position of this paper, in the quantum theory of liquids, is principally 
that of a critique of the contribution by Green. The relation of his 
formalism to wave-mechanics is discussed, and hydrodynamical equations 
are derived, in wavemechanical form, which are shown to be equivalent 
to Green’s results. A formula for the boundary stress is arrived at, and it 
is shown that the difference between the kinetic and thermodynamic 
pressures—on which Green’s explanation of the thermomechanical effects . 
in liquid helium IT is based—is identically zero, confirming the conclusions 
of Yvon and de Boer. 


$1. INTRODUCTION. 


THE system discussed here is that of N identical atoms (bosons or 
fermions), without effective internal degrees of freedom or spin, enclosed 
in a volume V and interacting by a potential U. H.S. Green (Born and 
Green 1947, Green 1948: referred to below as I and II respectively) has 
given the most complete account to date of the quantum mechanics of this 
model, and derived the hydrodynamical equations of motion in a suggestive 
form which demonstrates their correspondence to the classical equations. 
In this paper a treatment by conventional wavemechanical methods is 
given, and the results shown to be equivalent to Green’s. The virial 
theorem is derived, and the stress on the boundary identified in the surface 
term. A discussion of the thermodynamic pressure follows, in which 
de Boer’s (1949) conclusion is confirmed : that Green’s expression for the 
difference between the kinetic and thermodynamic pressures vanishes 
identically. 


§2. GREEN’s ForMALISM AND Wave-MECHANICS. 


We need a means of comparing results in Green’s formalism (I, § 2) 
with those derived below. For simplicity, consider a single particle only, 
with wave-function #(r) defined in a domain ® and satisfying a boundary 
condition B on its boundary, 8. If {¢,}, (v=1,2,...) is a complete’ 
orthonormal set of such functions, then 


Slt )= 2 Atretin) = eee 


* Communicated by the Author, 
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has the property 
ee 
| O(F,, Fo )f(r2)@r,—p(r,) in R, 


and obeys B on 8. 
In Dirac’s notation, if , is the wave-function corresponding to a ket 


|). 
p,(e') = <r’ |b). 


Hence by (1) crue a-0(r or’) (2) 
If A,= A(r,)=|1r1)<r, |, on re) 
then by (2), Aja 6(ts, fT), ERI R YT MENT > a ae (2) 
where r isa g-number. An operator A is expressed in the Schrodinger 
notation by 

«r’| PUL asta | meme. renew eh (A 
and in Green’s formalism by 

Gin oAs Wate) SAE Reus es obs 2.) say (6) 


Now if | £)<k |= p;, is the projection operator of a pure state h, 
pir a VP i 


(where the a, are real and $0) is von Neumann’s statistical operator 
(von Neumann 1943, Dirac 1947). By (5) and (6), if A and B are hermitian, 


[Ap,B\(r’, r’)= <r’ | Ap, B re 


=< (ki AAte aie 
=[Ad¢,(r ae ero ele at 
Hence E A,pB, | (r”, © j=0 (a) 
implies : eit Bie tS) 
co RA eat (0) | 


and conversely. 

All Green’s final c-number equations are in the form (8q@) (or an 
N-particle generalization, involving summations and integrations over the 
3N coordinates), and hence have equivalent wave-mechanical forms (8 6). 
For example, Green’s 

{p, P}= 3(pp+Pp) 
corresponds to Landau’s (1941) current operator 
Jin!) = 2(Ar)p + par), 
whose expectation 
Ce | Jee’) | k= BAL PvE PEP Walle 
is the conventional current at the point r’. Similarly 


Po) Te 
=) {{o, P} P} > — (WV PY wey) —C PPE Y NL. (9) 


mM 
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(The arrow denotes the correspondence (8), and the equality—for a pure 
state—(7)). Each of (8) is equivalent to 


Tr| pZAA(e) )B, |=o, 


and hence to the classical equation 


[[ pp. 2) 2Adp, 2)Bp, 2)8(e—©') dp de 


= { p(p. 8’) ZA,(p, r)BAp, FP p=0. . . . (10) 
In generalizing these results for N particles, r is replaced by the set 
f1,),...,¥%y. The effect of the symmetry restriction (Bose or Fermi) 


is that the p of (8) is in practice confined to the appropriate subspace of 
the space of {¢,}. The restriction is formally included by replacing all 
observables A by AQ, where Q=Q* or Q~ is the corresponding projection 
operator—which commutes with all A of physical interest. The corre- 
spondence (10) is not affected by the restriction. 


§38. THe Equations oF Morton. 
Write Schrédinger’s equation in the form 


mere, 
Hib(x,, Uy, ... , Xsy) = [—wA4+U (ay, Hg, ..., Xan) = th =; 


(11) 
where w=h?/2m, m is the mass of a particle, and 
A= = 0,0; 0; = 0/0x;. 

Let ka... . Xsy) = dih[bV p*—P* 7h]. + 84 One Rey 
Then 

V k= 20h, =(Jih (co) Hyp— pH] — — ml pry typed), 

| i 
Les 7 k=—m S (4). ae ES fee 
Also, : 


2 Fk =[(Hy)7 bP (Hy) (BY) AP ( 
= — 2p U— ol AYP p+ AY*T yp) — (UP (Ab*) +0 (4p) 
Hence, by the identity 
AUT P+ APT Y= 0 PWT PUPP (a. Pa, 
£ kU U=— of (Pup TY 


AV pt App Aye +27 yp . py} 
=— of [PPT P +O yl eA LAP). 


0 : vie Ts 
ak eV Ute (Ve Vp+V ov $*)—I7 7 (w*p)]=0. . (14) 
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§4. THe Macroscopic Equations. 
Rewrite the wave-function in terms of the coordinates of the particles : 


HOT age 2 ROS Rane mea 
where the r; are vectors in subspaces of the cartesian space of §3. If the 
surface S$ enclosing the particles is given by o(r)=0, then this 3N-dimen- 
sional space is bounded by the hypersurfaces S, : o(r;)=0. The boundary 
condition is taken as 

#=0 on the §,. ti iatiek ett 22.) (Loy 
We may also define a ** vector ” element of §;: dS, 

7S; 
(+0 that, for example, | r,.d8,=3VI1 Pr; ), 
j#t 

and vector operators V;, (V; .1;=36,;). The hamiltonian of (11) is then 


H=—w 2 4,4+U(rj, Yo, en eiens) ry). . . . . . (16) 


Let I,’ denote the operation 
hee = rer ere 
j#i 
and I; mean applying I,’ and dropping the suffix 7, so that 
I f(t)=VO f(r). 


etc., are defined similarly. Then 


ITV LV: ike) 
j 


ue 


aj? 


“Sj 
=1j(7,-kj)+ 21; | ky .48, 


JF 


: =1j(V,.)=V; .(1;'k)); by (15), 
Hence . Ties ki=sV2 (0K; ). 
Defining fee eee Lh) = (17) 
Ae . 
(13) becomes V eae Tm, Mt: OM Wr aan Way pe ore niet S) 


Since, apart from normalizing constants, J > mn, ui? (cf. Green’s I, 5.1) 
and n-> 7, (I, 2.7), (18) is equivalent to (I, 5.6) with g=1. Green’s 
more general cluster equations correspond to integration of (13) over 
fewer than (N—1) of the coordinates. 

Since (14) is a vector equation, we may write it 


@ 
Asi) | spect vu U | ae ats (Vy. Ty), .- + + + (19) 
where T,;=V; Vi b+V; PV > b*—3V,V ; (bp). 
Now, for j1, ie 
ieee T,)=1j| "dS, (0; °V b+V, WV, ¥*) 


since f=0, V;4=0, on 8;. 
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Defining 
R= aL V WV, p+ VW HI; 


S= lw S1,0,0; (bb) =40V Vn, 


(20) 
P= R—S, 
F= — 1, (p*47,U), 
(19) becomes I=F-7.P. 200 gis 2k Roti acted mee 


This result is equivalent to (i, 5.8) for g=1, since 
J > mn,u\”, F > mn,y\?, 
while for the remaining term of (I, 5.8) (there is a factor (1/m?) missing) : 
(1/m){{p1, pV }pT}> — fw ZT [WV Vip; V5 *)—VibVip+V pV *)] 
=v STV PV $+ Wi b*)—ViV, (bh) =P 


The equation corresponding to (I, 5.10) follows from (21) in the same way, 
given the definitions of v and d/dt, and the equation of continuity (18). 
The quantity P defined in (20) is the expectation, in Landau’s formalism, of 


()2 : (pi O(r, ri}, Pi 


which (by the considerations leading to (10)) reduces, in the classical limit, 
to the classical tensor of the momentum flux. Also S is proportional to 
the expectation of 


2[[p; d(r’, r,)], pi] (where [A, B]=(ih) (AB—BA)) 


oe 2 99 


and hence vanishes “‘ as h?”’, in the classical limit. 


§5. THe VirtaL THEOREM. 
From the definitions (20), the scalar of $ is 
S,= jw LT;A ,(p%)=gwodn, . . 2. . 2 . (22) 
and similarly 
R= FTA, (f"})— dab + pp*)]o dnp 2K, 
where 
K= Mier eent us eN Seek ere Rae Oe), 


is the kinetic energy expectation. Then from (21), 


a 
V (P.r)=(7 .P).r+P,=F.r—2J.r+2K+hodn. . (24) 
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vo 


For a stationary state, (0/dt)J=0. Then, on integrating (24) over the 
whole volume, ea 


s V 
[a.p.ra/ (RS Pdi)  . o eae (2b) 


where = —F . ris the virial of the potential U, and | : (An) d*r vanishes 


by (15). (25) is the quantum virial theorem for an energy eigenstate, or a 
statistical ensemble of stationary states. Putting 


a i<j 
the virial becomes 
PoAI,[r;. VV; Vilr) )b*b]+ 227 [ri (Valry) eh]. ©. (27) 
v U<J 
Or 
C=r.(VV,(r))n(r)+4fsVo(s)n.(r,r-+s)d’s, . . . (28) 
where 
NAL), Fg)=2ALT  (f*p). 
V<) 
By (15), on the boundary 
V p=n,0,;4, 
where n; is the unit (outward) normal to §8,;; and hence 
onS, P=nnw2I,(0,;$%0,;$). . . . . - - + (29) 


The interpretation of (29) as the tensor of the stress on the boundary §, 


suggested by the form of (25), is supported by the following considerations 
for one dimension. 


The Schrédinger equation in one dimension is 


| 0 Fa + U(w)—E | bla) —0, ee Pr 0) 
where 
Lim U >E. 
@>+ 0 


Suppose that, for c>0, the potential has the form 


LEA, U=0, l | ie 
x>Aa, U=E,=constant, | 
where E,>E. Then the solution of (30) for 7>0 is 
a<a, p=Asin [k,(x—a)+4], 
“>a, YA sin (¢) exp [—k,(a—a)], | (32) 
wk=KH, wk2=EK,—E. | 


Since (d/da)y is continuous, 


E ae 
tan (¢)= — nlte=— (qa) iia . f 5 (33) 
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The contribution of the “ potential step ” (31) to the reaction of the system 


1s 


B= Lim |" py dU=EWW 0 


£0 
ehek ad \ (a ; 
=| AE, sint (6) =00(7. 0). (aaa) (lo tan? (4)/B), 
If E is fixed while EK) > ©, 
(a) s>0, nS a, 


and hence 


which is the one-dimensional analogue of (29). 


§ 6. APPLICATION OF FROHLICH’S PERTURBATION THEORY. 


The interpretation of (29) suggested above, is confirmed by the theory 
of boundary perturbations proposed by Frohlich (1938). In the notation 
of (11), let % satisfy HA=KHys in ® and vanish on its boundary §, and let 
ob’ satisfy Hb’=E’p’ in ®’ ]@ ® and vanish on 8’. Then by Green’s 
theorem 


(H’—B) | Wylde =—w | Y*d,p' dS. S36) 

Choose 7’ so that, as 8S’ > 8, %’ > % and E’ > E. Suppose that, if x is on 
8, x’=x-+Az(x) is on 8’. Then it is assumed (Frohlich 1938) that 

Anab'(x+ Az) =Oyah(x)+O(X) -. 2 2. . . (89) 


(i. e. # “ moves with the boundary ”’, to the first order in A). Then (36) 
gives ; 


W=B—Aw { B,W\(2,h)2AS+OO2). 188) 


In the notation of §4, if the surface S enclosing the particles is perturbed to 
S’ by r’=r--As(r), then (38) integrates to 


dik ‘S ‘ 
Car)... =~o| | 21, nit Ani) | s.dS=— | dS SPAS: , , (39) 


1 Ct tw £ j a = fe » ay 1 A nf “ 
Conservation of energy then requires P to be the tensor of the stress on 
the boundary in a slow perturbation, provided Ehrenfest’s theorem applies. 
For a Gibbs ensemble, the free energy F is given by 


fee) 
eta = et. Daren stares Sie or. NE) 
y=1 
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where E, are the energy levels of the eigen-functions of H, of Bose and 
Fermi symmetry, and 8 1=kT. Then 


oF loge = di v — BR 
i), =Z, i io va =—f ds. 2-420 | $847 ¢...2 6) 
- and therefore p= a, 


the average of P over the eR distribution, is the boundary stress 
for isothermal displacements. For a simple dilation 


yo _(t\4 
s=rf, aVee 3 Ate 6. ag Eo et BAe (42) 


Substituting (42) in (25), (39) and (41), the thermodynamic pressure is 


dk. 
p=- = (Gy) iy [zrzemeK,-1) | der. 


-V = ane. 
Ie. BpV = | (SRST ior tay ty ee Cy sy ee (43) 
There is an arbitrariness in (39), which can be written 
dK 
— Fx | 48. @R—AS).s 


where 24—fB=1. For the case (42), (43), s=r we then find 


“V a V 
oaks | eK Tyres (7 *) [ Side, a0. 


~ an 
“V. 
Since S.d°*r=0 by (15), we then retrieve (43). 


Comparison of (23), (28) and (43) with (I, 3.21) and (I, 3.22) shows 
that p is numerically equal to Green’s p,, the first term in his series for p. 
Thus the quantity 7=p,—p, which Green identifies (II) with the energy 
associated with the thermomechanical effect in liquid helium II, appears 
to be zero. However, this conclusion has been the subject of a recent 
controversy (Yvon 1948), (Massignon 1949), (de Boer 1949), (Green 1949), 
(Horowitz 1950). A further discussion of the question is therefore 
attempted in the following paragraphs. 


§ 7. APPLICATION OF BRILLOUIN’S THEORY. 

The result (39) may also be obtained by the independent perturbation 
method originated by Brillouin (1937). For simplicity, the work is given 
in terms of the case N=1. The particle is enclosed in the surface 8, 
on which the wave-function (r) vanishes. Let S be continuously 
deformed to 8S’ by r’=r-+As(r), and choose s to be normal to $. Then 
s(r) may be chosen inside 8 so as to agree with the boundary value, be 
continuous, and satisfy 

Vees—constant, s=Vdi. .°. . .». (44) 
3U2 


956 P. J. Price on the 


The choice is unique, and has the property 


Assos OS eee 
If H'b(r)= Heb (r’) = Hy (r-+As), 
then H’=H-+Af{w[2 ANN ae 5. 
Or, by (45), 
dH. 
a= (5) =2u(V.s):V 0 +[s.V Ut ores GO) 
dX}, 


Since, by (44), the eigenfunctions, ,, of H remain orthogonal under the 
deformation, the change in the eigenvalues is given by 


dh, : 
Ty = (be LG). 1 aetna Ew eee a 
Now : 
V.(P.s)=(V.P).s+P:Vs 
= —p}(VU) .s +o’ s : [Vee b7 ot IVP (bp) 
for a stationary state. Then, integrating by parts and using (45) and the 
boundary condition, we find 


78 V V 
—| sags Se [V s:b*V Pp]der 


ae 
WIGIW=G) lH... 2. 2... (48) 


The generalization for N particles is om straightforward. Thus (39) 
is confirmed, in a form which covers degenerate unperturbed eigenstates. 
If the boundary wall exerts a short-range force on the molecules, with 
potential U, which “‘ moves with the wall” in the perturbation, then the 
derivation of (46) shows that the boundary stress is still given by (46) 
and (47) if s.7U is replaced by 


s.V(U,+U, )4e(Se ans, a eed 


where n is the unit normal to 8, and —zn is the vector from rto S. For 
the dilation case, s=r, (49) reduces to r.VU, and (43) is still valid. 


§8, THe THERMODYNAMIC-PRESSURE SERIES. 


Derivations of (43) have been given (de Boer 1949), (Green 1949) which 
make use of the transformation of the partition sum Z of (40) to 


#—Lim(p, n) 2 (ko, |) . yg UN OO 
where Git « A ly spook arenes poate eee EN CLO] 


{| &)} is any complete (in the Tans or Fermi subspace) normalized 
orthogonal basis, and Lim(p, ) means 


Lim Lim 
p> oo n—> 0 
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taken in the order from right to left (it is, of course, not permissible to 
interchange the limits in (50)). By (46), 


d : r 
7H =HG+HOGHS ... +GH= 2 (eee Geo, (52) 


es 
where GHG, ,—G,.,H, °G,=G: 
Hence, by (50), 
Peete hes | caus (kh) GyHe-[ ke) |. (63) 
dr k=1r=1t=1 t!(r—t)! )! a 4 


If we choose, for the | k), the eigenkets of the unperturbed hamiltonian, 
then 
(k| GH" | k)=<k| HG, | EY —(k| G_, | &)E, 0. 
If, H being hermitian, we put 
(k| HG, 1 | k)=E,¢k | G1 | *) 
(it is pointed out below that this step requires a special justification), 
then Ge Ga yea On eb eels cc 1 - ca MOL) 
Substituting (54) in (53), by (48) we retrieve (41) and (43). 
Since (for s=r) 


G3 V = (i) (Hp ar—p.rH), .. .« «-.. (5b) 


the argument leading to (54) also requires (k| GH"|k) to vanish. It has 
been pointed out, in this connection, by Green (1949) that expressions 
of the form 


n v 
THAB—BA)=Lim(n) 2 2 (A,,B,,—B, Aor 

r=loe=1 
need not vanish, with a dense or continuous basis {|7)}, if (r|AB|7r) 
is not uniformly convergent with respect to 7. There is, however, another 
approach to such questions, which was pointed out by Fuchs (1940). 
H is specified not simply by a differential operator, such as (11), but with 
respect to the linear space of a basis {| k)}—+. e. to the domain symmetry 
and boundary condition of the corresponding wave-functions. H_ is 
then only hermitian in virtue of this boundary condition, and if A |) 
does not satisfy the latter (A is then called by Fuchs “ inadmissible ’’) 
we cannot assume that 


(J | (HA | k))=[¢h | (AH|)>)1* 


Differential relations such as (46) and (55) are still correct, but deductions 
by matrix algebra may be invalid. For example, a virial theorem may 
be deduced as follows : 


, (¥; + Pi Ps - Ti) 
==") Ee eee eiman FU, 
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and hence for a stationary state 

«| Si(p?/m—r,; .V,U)| >=9. 
The surface-integral term is missing because, although the commutator 
relation [1,, p;]=1,5,; is applicable to r; and p; as differential operators, 
the Heisenberg equation 


d 
ShIAl =A: HI) 


is then invalid where A is not admissible. Now, with the boundary 
condition s=0 which we are considering, G and the G, are admissible, 
since G—2H reduces on the boundary to a function of x, but p. x= seap 
isnot. It follows that the diagonal elements, in the energy representation, 
of G,, G,, . . . necessarily vanish, but not those of G. 


§ 9. THE COORDINATE REPRESENTATION, 


Green’s (1949) treatment is based on the Schrédinger representation 
of (50) (Slater 1931, 1933): 
eM ian ee ed hae a etm eng TS 
A typical term in the resulting series for dF/dA (simplified by not 
employing the relation (55), and not symmetrizing the operators in each 
term) is proportional to 


| d2Nx'[G,eP2](x’, x’)= | dx! 3 fb, *(x')[Ga,(x le}, - (57) 
v=1 


where the #, are the eigenfunctions of H. For t1, since | %,| is bounded 
and (at least in the dilation case (42)) G, reduces to a finite function of 
position, (57) is 


Er eB, | [Oxf] | = Ee? Qy| Gy] vy 0. 
v=1 v=1 

Thus de Boer’s (1949) statement, that Green’s series reduces to its first 
term, is confirmed. A recent paper by Horowitz (1950) appears to 
suggest : 


« 


(a) that G is inadmissible, in Fuchs’ sense, 
(6) that the reasoning, corresponding to (53), which leads to Green’s 
series. is therefore invalid. 

We have seen that G (though not Green’s operator) is admissible for 
the boundary condition 4=0, which Horowitz also adopts in the context 
of his criticism (G is not admissible for a periodic boundary condition, 
which leads to difficulties if applied in this problem). But in any case, 
from the present point of view, the admissibility of G affects, not the 
validity of the derivation of Green’s series, but the vanishing ov non- 
vanishing of the terms of this series. However, the validity of the use 
of (50), in the form (56), needs examination—since the transformation 
matrix 


Ty, = «x | vp) 
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is singular (T'T(x’, x”)=8(x’, x”). 


iz I<, n, Re | dx! 5 y oad *(x’)[H", (x’)] 


and Lim (S) means the limit G + ¥, then of the six expressions 


Z,=Lim (n, 6, d)I, ] 

Z,=Lim (n, A, S)I, 

Z,—=Lim (A, n, S)I, | 

Z,—Lim (A, S, n)I, me 
Z;,=Lim (6, A, n)I, 


Z,= him’ (S,'n, A)1; 


Z,, Z, and Z, are divergent. Using the fact that the |1,| are bounded, 
f= /4=7,,.. But 


Pee eg eR ee oe (61) 
yv=1 
while Z, is evidently the Slater (1933) sum, (56). Then, by (52) and (59), 
dZ_ dZ, 
dX. dx 


pate x')[G,_ HOY, (x 


Lim (2, A, n) | Ox E13 z en a ee x’)[G,_ 1b, (x ee 
2 y=1lr=1 t=1 


=lrm-(©, 1; ny |[~ BNx! 3 3 pallial 


PO et le 


Wy 


(62) 
Since 
2 < eae Hsp" (x) (Gy, (* 
s=0 98 ! 
converges uniformly with ce: to t, x’ and y, 
Jugueret i 
aa e- wah iby (X ‘LG, 1h, (X x’) ao See (63) 


which also converges ne with respect to x’ and v, and hence 


Lim (n){ rte 


dZ coe t 
—~ = Lim(<, r) | deux’ z aa aad (Gren! (x 2.) head Oe) 
dr es 
Interchange of the remaining two limits results in 
dL Se(—P) (Bn —pH : 
=, = (— ———_— Lex’ (Gye Xa as OD 
TA (—f ot 71 | d a i NGxeo Ka (65) 


which is essentially Green’s series, and which we have shown to reduce 
to its first term, 
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CONCLUSION. 


Of the mathematical results in (I), the work in this paper confirms the 
dynamical equations (5.6), (5.11), and thus vindicates the formalism 
which led to them, but is at variance with the more controversial 
conclusion of §3: that the thermodynamic stress on the boundary 
(or on a “probe” inserted in the liquid) differs from the kinetic 
stress. However, it is not clear how far the theory in (II) really 
depends on (I). The formuls for the thermodynamic stress necessarily 
refer to the whole volume of the system, the integrals over which involve 
a boundary contribution associated with S, (20). The theory given in (IT) 
of the thermo-mechanical effects (cf. II, 3.20) must be regarded as an 
empirical one (7 was in fact computed from Kapitza’s experimental 
results), while the discussion of the normal viscosity (I, 6.1) (II, §2) 
seems to depend on a conception, of the local order in the liquid, 
whose foundation is not clear. Further progress will probably depend on 
the development of a theory of the structure of the quantum liquid when 
almost in equilibrium. It may be, for example, that there is a long-range 
order which makes it meaningless to speak of local quasi-isothermal 
processes, to which results similar to (41) would apply. 


I wish to express my thanks to Mr. H. N. V. Temperley, of the Royal 
Society Mond Laboratory, for his advice and encouragement, and I am 
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LXXXVITI. CORRESPONDENCE. 


Latitude Effect of Cosmic Ray Stars. 


By 8. Larrrmorg, 
Imperial College of Science and Technology, London *, 


[Received July 25, 1950.] 


Inrorp 200 G.5 emulsions were exposed vertically for 27 days at 
altitudes of 11,550 ft. and 41,400 ft. at a latitude of approximately 12° S., 
in the Bolivian Andes. The plates were exposed under a corrugated iron 
roof. Since the time required for transporting the plates to and from 
Bolivia was considerable, control plates of a similar type to those actually 
exposed were used to find the number of stars collected before and after 
the exposures. . 

The results are given in Table I., together with figures for the 
Jungfraujoch. The latter were calculated from results at 11,500 ft., 
obtained using 200 ~ G.5 plates, assuming a range of 145-10 gm./cm.? 
of air for the star producing radiation. 


TABLE I. 
No. of stars (< 3 prong)/c.c./day. 
Height. Ratio 
(feet) Eades pone Jungfrau/Andes 
11,550 TT7+41°5 16-4-L0-5 2-1 +0-4 
14,400 13-4+1-7 27-5+1:3 2-05 -+0-27 


The value of 2-1+0-4 for the ratio of the number of stars with three 
or more prongs at the Jungfraujoch to the number at 12°S. is in fair 
agreement with results of Beets, Morand and Winard (1949) who found 
a value of 2-4 for the ratio between 44° N. and 1°S. at a height of 14,700 ft. 

In Table II. are shown the size distributions for all stars at the Andes, 
and those at the Jungfraujoch, normalized so that the ratio of all stars 
with three or more prongs at the Jungfraujoch to all stars with three or 
more prongs at the Andes is equal to 2:1. As can be seen, the latitude 
effect decreases as the star size increases. 


TABLE [I]. 
No. of Stars. Lead 5 
: Ratio. rau 
Star size Andes Jungfraujoch ato rerne 
3,4 193 +25 498 +22 2-6 +0-35 
D6 83:5+14 161-13 1-93--0-35 
Tn tes 42 19 52+ 7 1-24-++0-32 
~ Total 360 755 2-10 


* Communicated by Professor Sir George Thomson, F.B.5, 
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The figures in Table II. are plotted graphically in fig. 1 and are shown 
to lie on a smooth curve. 

Now the magnetic cut-off energy at the Andes is 15,000 MeV. and at 
the Jungfraujoch is 2000 MeV., and it is possible to estimate from the 
ratio of the number of stars at the Jungfraujoch to the number at the 
Andes the proportion of the stars at the J ungfraujoch that are produced 
by primaries with energies between 2000 and 15,000 MeV. 


Fig. I: 


4 5 6 7 8 9 10 i I2 IS 


Star size. 


Let x be this fraction. Then 


2-l=1/(1—2) oy Motroee 


a? 


i.@. only one-half of the stars at the Jungfraujoch are produced by 
primaries with energies between 2000 and 15,000 MeV., although, 
assuming an K~* energy spectrum for the primaries, there are approximately 
50 times as many primaries in this range as there are with energy greater 
than 15,000 MeV. This result is, at first sight, rather surprising, since 
most of the stars have energies below 300 MeV. However, the majority 
of the stars are produced by neutral particles, probably neutrons, which 
must have been produced by the primaries high in the atmosphere, and 
further, the work of Harding (1949) shows that a star producing particle 
can probably produce more than one star. It would, therefore, seem 
probable that a primary produces several particles high in the atmosphere 
each of which can produce several stars, in which case the lower energy 
primaries would become less important, 


Correspondence 963 


ACKNOWLEDGMENTS, 


We have pleasure in expressing our gratitude to Mr. Clifford Waite of 
Consolidated Tin Smelters Ltd., and to Mr. D. C. Deringer, Jr. of Platino 
Mines and Enterprises Cons. Inc. for exposing the Andes plates used in 
this experiment. Our thanks are also due to Sir Roger Makins, K.C.B., 
of the Foreign Office, who arranged for the transport of the plates to and 
from Bolivia. We are grateful to Professor Sir George Thomson for 
interesting discussions on this work. 


Imperial College London. S. LATTIMORE. 
REFERENCES. 
Beets, C., Moranp, M., and Winarp, L., 1949, C.R. Acad. Sci., Paris, 229, 
1227-8. 


Haropina, J. B., 1949, Phil. Mag. (7) 40, 530. 


Note on the Heavy Nuclei of the Cosmic Radiation. 


By A. D. Darnton and D. W. Kent, 
H. H. Wills Physical Laboratory, University of Bristol *. 


[Received July 30, 1950.] 
[Plates XXX. and XXXI.] 


It is now well known that the atomic nuclei among the primary cosmic 
radiation are rapidly absorbed in the atmosphere as a result of nuclear 
interactions and loss of energy through ionization (Bradt and Peters 
1948, 1950, Freier ef al. 1948). Further, in the nuclear collisions which 
they make, these primary particles are often not absorbed catastrophically. 
They frequently suffer fragmentation so that fast heavy particles of smaller 
mass emerge from the encounter—see Pls. XXX. and XXXI._ It follows 
that the distribution in mass and energy of the heavy particles, as well 
as their intensity, changes rapidly as they penetrate towards the earth. 
Recent experiments in this laboratory with plates exposed at 68,000 feet, 
showed the presence of substantial numbers of heavy particles at this 
altitude. The observations have allowed us to determine the character- 
istics of the flux of heavy particles at this altitude. The deduced values 
of the intensity of the heavy nuclei at the top of the atmosphere, and their 
mean free path, are in good agreement with the results of other workers. 
The exposure was made with 400 Ilford G5 emulsions by means of 
free balloons at a geomagnetic latitude of 54°.N. Level flight was 
maintained at 68,000 feet (52-+1 gm./em.2) for 100 minutes. Since very 
few heavy nuclei are observed at altitudes lower than 55,000 feet, and since 


i r BAS 
* Communicated by Professor C. F, Powell, F.R.S, 
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the time of ascent and descent above this altitude was small compared 
with the time of level flight, it was a sufficient approximation to assume 
that the observed nuclei entered the plates at 68,000 feet. 

Our observations were made on 118 nuclei with charge Z>5. The 
charge of each nucleus was determined by the delta-ray method described 
by Bradt and Peters (1948). The number, v, of delta-rays having an 
energy greater than W,, produced per unit length of its trajectory by a 
particle of charge Z, is given by: 


where fc is the velocity of the incident particle. Using the range-energy 
relation, the charge can thus be determined within an error of plus or 


Fig. 1. 


Number of nuclei. 


tein waite 


0 
Charge Z 


Distribution in charge of 118 nuclei. 


minus one unit in the majority of cases, if the length of the track is suffi- 
ciently great. In cases where particles of low charge cannot be observed 
over a long range, the uncertainty in the charge determination is greater, 
The spectrum determined in this investigation is shown in fig. 1. 

In order to obtain a value of the flux at the top of the atmosphere it was 
convenient to consider only those tracks which passed through the 
emulsion of the outside plates of the stack. The angle of each track with 
respect to the zenith was determined from the dip and the angle between 
the projection of the track on the plane of the emulsion and the vertical 
side of the plate. The observed distribution was then corrected for 
absorption in the stack and for geometrical orientation of the plates. 
The number of tracks observed having a zenith angle between @ and @+-dé 
is given by 


dn=I(6).t.A.4sin20. (6). dd, 
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where ¢ is the time of exposure, A is the area surveyed, and 4 sin? @ . f(0) 
is the correction factor due to plate orientation and to the choice of 
minimum track length (500 in this investigation). The function 


f(@=1— | 1— (<3) | 


where d is the depth of the emulsion and / is the minimum track length 
considered, is operative only when the zenith angle exceeds sin-! d//=53°. 
Assuming the absorption in the atmosphere to be exponential, 1(0), the 
intensity of the radiation at an angle 0, is related to the flux at the top of 
the atmosphere by the equation 


I(0)=1, exp [(L/A) sec 6], 


Fig. 2. 


0 
Secant 6 
Angular distribution function. 
L being the distance in gm./cm.? from the level of the exposure to the top 


of the atmosphere, and A the mean free path for disappearance in air. 
We can express equation (2) in the form (Hoang Tchang-Fong 1950) 


log, | 5/s ar) | loge (dey. A bh sec 6. 
A plot of this function, shown in fig. 2, gives a value —L/A for the slope and 
- for the intercept at sec 0=0, log, (4.1).A.¢). From these two values we 
obtain for the mean free path 
A~28 gm./em.”, 


and for the flux Ij~1:3 x 10°° particles cm. sec.~! steradian™?. 
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It is probable that an appreciable fraction of the nuclei observed at this 
altitude are not of primary origin. Of eight observed collisions caused by 
nuclei of charge Z> 6, four cases resulted in a reduction in the charge of 
the primary, but little change in its velocity and direction. 

Tn 34 cases in which an increase in the delta-ray density with range was 
observed, it was possible to determine the kinetic energy of the nucleus. 
The values thus obtained were corrected for ionization losses in the air to 
determine the corresponding values of the kinetic energy at the top of 
the atmosphere. Since collisions involving degradation of the charge 
occur, it is reasonable to assume that some of the particles are splinters of 
larger, more heavily ionizing nuclei. The values obtained, which lie 
between 0:45 BeV. and 2:0 BeV. per nucleon, therefore represent lower 
limits for the energies of the particles at the top of the atmosphere. 

The present experiments are being continued with the application of 
new techniques. 


The authors wish to express their gratitude to Professor C. F. Powell, 
F.R.S., for making this investigation possible, and to Dr. C. Franzinetti 
for helpful guidance and discussion. 
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Loss of Energy of Fast Protons in Matter with Particular Reference to their 
Range in Carbon. 


By ‘A. E. Taytor, : 
Atomic Energy Research Establishment, Harwell, Berks.* 


[Received August 3, 1950.] 


In using the 110-inch cyclotron at Harwell for measurements of total 
scattering cross sections and neutron energy spectra (Taylor, Pickavance, 
Cassels and Randle 1950) carbon absorbers of varying thickness have been 
used in a triple coincidence telescope detecting protons scattered from a 
polythene radiator. The range in carbon was used to deduce the energy 
of the protons, since carbon was a convenient material having a low atomic 
number. This choice was made in order to minimize multiple coulomb 
scattering and nuclear absorption for efficient operation of the detecting 
telescope. It was therefore required to compute the range coy relati 

in ane for protons up to 260 MeV. Be ae es eae 
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A nucleus of charge Z=13-+1 with a velocity 8B=0-55 suffers a collision in the 
emulsion and gives rise to a “‘jet’’ of six a-particles (considered to be com- 


ponents of the primary particle) and low energy disintegration fragments of 
Observer: Miss J. HeRapErR, 


the target nucleus- 
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A fast nucleus of charge Z—16-+-1 interacts with one of silver or bromine in the 
emulsion. The nucleus emerging from the encounter has a charge Z=9-- 1 
and is accompanied by twenty-five singly charged particles moving at 
relativistic velocities, and in directions contained by a cone of semi-vertical 
angle 30°. The tracks of fifteen of these particles have been reproduced in 
this photomicrograph. Observer: Miss EB. JAMEs. 
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The calculations have been made using the theoretical equation for 
protons given by Livingston and Bethe (1937) 


dK 47Net[_, 2mv? z 
tame |! py PCs | 


where —dE/da is the rate of loss of energy, Z’ the effective number of 
electrons per atom, I the mean excitation energy of the atom of the 
absorber and Cx a correction term. 

For elements of low atomic number, Z’ is taken as the total number of 
electrons, Z, and the non-excitation of the inner shell electrons for low 
proton energies is corrected by the term Cy. Values of Cx are given by 
Livingston and Bethe (1937) in terms of the energy of the proton and a 
critical energy depending on the ionization potential of the electrons. 
For elements of high atomic number Z’ is taken as the effective number of 
electrons (Hénl 1933). 

The value of I is usually chosen to give agreement with experimental 
values for the energy loss! at low energies. On the basis of the Fermi 
Thomas model of the atom, Bloch (1933) showed that approximately a 
linear relationship should hold between I and Z for elements of high 
atomic number. The constant of proportionality, «, given by Wheller 
(1941), is 11-5 eV. while Mano (1934 a, b) has given the value of 10-5 eV. 
for elements with Z greater than 20 but shows that « fluctuates for Z less 
than 20. The values adopted by Livingston and Bethe (1937) for hydrogen, 
helium, aluminium, air and gold differ slightly from those of Mano, but 
the variation of « with Z follows the same pattern. 

Before the computation of the range energy relation in carbon can be 
carried out, an appropriate value of I must be found. In view of the 
variation of « for elements of low atomic weight, it was decided to find a 

‘value which gave an atomic stopping power relative to aluminium in 
agreement with experimental values, rather than take the theoretical value 
of 11-5eV. Using for Z’ the total number of electrons both for aluminium 
and earbon and a value of I for aluminium as 150 eV., the value of I for 
carbon was found to be 50-3 eV., and hence « is 8:4 eV. This value seems 
to fit in with the pattern of variation of « with Z. Similarly for silver and 
copper the calculated values of I were 365 eV. for silver and 293 eV. for 
copper, using values of Z’ 39-6 and 27-4 respectively. 

Using the value of I as 50-3 eV. for carbon and the appropriate values of 
C,, dE/dx was calculated at 20 MeV. intervals up to 200 MeV. Below 
20 MeV. the known experimental values were taken. The range was then 
computed by numerical integration and is given in Table I. 

These values for the range have been used as a measure of the energy of 
protons scattered into a triple coincidence telescope from a polythene 
radiator bombarded by neutrons ejected from beryllium and carbon 
targets. The end point of the spectrum of the scattered protons was 
compared with that using aluminium for the absorber. ‘The energy in 
this case was deduced from the range energy curves of Smith (1947). 
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TaBeE [. 
1 dB | 
Energy p dx | Range 
MeV. MeV./gm. cm.>?|_ gm. cm.~? 
25 132-5 | 0-013 
5 17-3 0-0374 
10 44-15 0-127 
15 31-58 0-264 
20 24-96 0-442 
40) 14-16 1-55 
60 10-23 3°24 
80 8:16 5-45 
100 6-89 8-13 
120 6-03 - 11-24 
140 5:39 14-76 
160 4-91 18-65 
180 4-54 22-89 
200 4-23 27-45 


The energies so determined agree with each other to 2 MeV. and to within 
2 MeV. of the value expected from the theoretical energy of the neutrons, 
the experimental error being also about 2 MeV. 


August 2nd, 1950. A. E: TAYLor. 
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